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Abstract 


Experimental  and  theoretical  studies  of  hypersonic  boundary-layer  stabilization  using  a  passive 
ultrasonically  absorptive  coating  (UAC)  of  regular  microstructure  were  conducted  in  order  to 
evaluate  the  UAC  performance  for  hypersonic  laminar-flow  control.  The  coating  is  a  porous 
sheet  perforated  with  equally  spaced  blind  cylindrical  holes  of  the  average  pore  diameter  50  pm, 
depth  450  pm  and  average  spacing  100  pm.  This  coating  is  similar  to  that  tested  in  the  GALCIT 
T-5  shock  tunnel  of  California  Institute  of  Technology.  Linear  stability  problem  for  hypersonic 
boundary  layers  on  the  porous  wall  was  formulated  using  the  boundary  conditions,  which 
simulate  absorption  of  disturbance  energy  by  a  porous  layer.  These  boundary  conditions  result 
from  modeling  of  acoustic  disturbances  propagating  inside  a  long  cylindrical  pore.  Solutions  for 
the  dynamic  density  and  dynamic  compressibility  are  refined  to  account  for  the  gas  rarefaction 
effect  associated  with  Knudsen  layers  on  the  pore  walls.  Parametric  calculations  of  the  UAC 
effect  on  the  second-mode  stability  are  performed  for  a  sharp  cone  at  freestream  Mach  number 
=  6 .  It  is  shown  that  the  UAC  strongly  stabilizes  the  second  mode  in  the  cold-wall  case  typical 
for  practical  applications.  Benchmark  measurements  of  acoustic  absorption  by  the  porous  coating 
are  conducted  using  the  acoustic  interferometer  developed  in  ITAM.  These  measurements  agree 
well  with  theoretical  predictions  of  the  UAC  impedance  and  absorption  coefficient.  The 
theoretical  model  gives  a  similarity  law,  which  allows  for  extrapolation  of  benchmark  data  to  low 
pressures  and  hi^  frequencies  typical  for  experiments  in  hypersonic  wind  tunnels  and  flight. 
Stability  experiments  are  performed  on  a  7“  half-angle  sharp  cone  in  the  T-326  wind  tunnel  at 
freestream  Mach  number  5.95.  The  cone  has  a  longitudinal  half  of  its  surface  solid  and  the  other 
the  porous  coating.  Hot-wire  measurements  of  natural  and  artificially  excited  disturbances  are 
conducted  on  both  porous  and  solid  sides  of  the  model.  Analysis  of  mean  profiles  and  RMS 
mass-flow  pulsations  show  that  the  UAC  roughness  does  not  cause  premature  tripping  of  the 
boundary  layer.  Natural  disturbance  spectra  indicate  that  the  second  mode  is  a  dominant 
instability.  The  UAC  stabilizes  the  second  mode  and  weakly  affects  the  first  mode  that  is 
consistent  with  the  theoretical  predictions.  Measurements  of  artificially  excited  wave  packets 
show  that  the  UAC  leads  to  substantial  decreasing  of  the  wave-packet  growth.  The  experimental 
data  on  phase  speeds  and  amplitudes  of  the  second-mode  disturbances  are  compared  with 
theoretical  predictions.  Satisfactory  agreement  is  obtained  on  both  solid  and  porous  surfaces.  The 
linear  stability  theory,  benchmark  measurements  and  wind-tunnel  experiments  confirm  the 
concept  of  hypersonic  boundary-layer  stabilization  using  passive  porous  coatings  of  regular 
microstructure.  Theoretical  and  experimental  results  of  this  effort  can  be  used  for  evaluation  of 
the  UAC  performance  and  design  of  hypersonic  laminar  flow  control  systems. 
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Nomenclature 


For  Chapters  I  and  III 

A  =  disturbance  amplitude 

A^,A^,  A^  =  admittance  components 

A^  =  thermal  admittance 

C  =  dynamic  compressibility 

C'p  ~  specific  heat  capacity  under  constant  pressure 

Cx  =  longitudinal  phase  velocity  of  the  wave  packet 

/  =  frequency 

F  =  a;  V* frequency  parameter 

h  =  porous  layer  thickness 

Kn  =  Knudsen  number 

r  /ul 

M  =  Mach  number 

m  =  propagation  constant 

p  =  pressure  disturbance 

P  =  mean  pressure 

Po  =  stagnation  pressure 

Pr  =  Prandtl  number 

Q  =  mass  flow  disturbance 

fp  =  pore  radius 

R  I  vl 

Re  =  Reynolds  number 

Rei  =  unit  Reynolds  number 

s  =  pore  spacing 

SA,  SF  =  amplitude  and  phase  spectra  with  respect  to  the  transversal  wave  number  P 

SAo  =  value  of  SA  at  P=0 

T  =  temperature 

To  =  stagnation  temperature 

u,v,w  =  velocity  disturbance 

U,V  =  mean-flow  velocity  components 

X,  y,  z  =  Cartesian  coordinates 

X,Y,Z  =  hot-wire  coordinates 

=  characteristic  impedance 

t 

a,l3  =  wavenumber  components 

6  =  boundary-layer  thickness 

S*  =  boundary-layer  displacement  thickness 

£  =  T  /  n  small  parameter 

7  =  specific  heat  ratio 

p  =  density 

p  =  dynamic  density 

K  =  heat  conductivity 
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For  Chapter  II 


C 

c 

d 

f 

H 

h 

k 

K 

k, 

L 

n 

P 

Pr 

q 

fo 


s 

h 

W 

z 


=  complex  compressibility 
=  speed  of  sound 

=  diameter  of  the  interferometer  tube 
=  frequency 

=  r^^p^u/rj  scale  parameter 

=  porous-layer  thickness 
=  wave  number 
=  structural  constant 

=  added  mass  coefficient 
=  length  of  the  interferometer  tube 
=  standing  wave  faetor 
=  pressure 
=  Prandtl  number 
=  pore  length 
=  pore  radius 

=  thickness  of  the  porous  sample 
=  shape  factor 

=  wave  resistance 
=  specific  impedance 
=  input  impedance 


a  =  absorption  coeffieient 

p  =  complex  dynamic  density 

CO  =  angular  frequency 

^  =  porosity 

<j  =  flow  resistance 

rj  =  vicsosity 

cp  =  absorption  phase 

X  -  wavelength 

Y  =  specific  heat  ratio 

r  =  propagation  constant 


List  of  Acronyms 

ADC  analog-digital  converter 

CCA  constant  current  anemometer 

ITAM  Institute  of  Theoretical  and  Applied  Meehanics,  Russia 

PSE  parabolized  stability  equations 

RMS  root-mean-square 

UAC  ultrasonically  absorptive  coating 
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Introduction 


Laminar-turbulent  transition  leads  to  substantial  increase  of  the  aerodynamic  drag  and  surface 
heating  of  hypersonic  vehicles  (Lin  et  al.,  1984;  Tartabini  et  al.,  2000).  Transition  also  affects 
engine  performance  and  aerodynamic  control.  The  ability  to  increase  the  laminar  run  is  of  critical 
importance  in  design  and  optimization  of  aerospace  planes.  Estimates  for  the  National  Aerospace 
Plane  (Whitehead,  1989)  indicated  that  payload-to-gross-weight  ratio  would  nearly  double  if  the 
boundary  layer  on  the  vehicle  surface  were  fully  laminar,  compared  to  fully  turbulent. 

Aspects  of  hypersonic  boundary-layer  transition  control  are  discussed  by  Kimmel  (2003). 
Smoothing  and  shaping  of  the  vehicle  surface  help  to  avoid  early  transition  due  to  roughness, 
leading-edge  contamination  as  well  as  cross-flow  and  G"  rtler  instabilities.  However,  with  these 
measures  the  laminar  run  may  still  be  short  because  of  the  first  and/or  second  mode  amplification 
(Mack,  1984;  Malik  et  al.,  1990).  The  wall  cooling,  which  naturally  occurs  on  hypersonic  vehicle 
surfaces,  strongly  stabilizes  the  first  mode  (Lysenko  &  Maslov,  1984)  while  destabilizes  the 
second  mode  (Malik,  1989).  This  indicates  that  hypersonic  laminar  flow  control  concepts  should 
address  the  second-mode  instability,  especially  for  cruise  airbreathers,  due  to  the  extended  period 
of  hypersonic  flight  at  relatively  small  angle  of  attack.  Since  these  vehicles  have  predominantly 
planar  shape  with  sharp  leading  edges  and  the  local  Mach  numbers  are  high  (M>6),  the  second 
mode  may  be  a  dominant  instability  leading  to  laminar-turbulent  transition. 

The  second  mode  is  the  result  of  an  inviscid  instability  present  due  to  a  region  of  supersonic 
mean  flow  relative  to  the  disturbance  phase  velocity.  This  mode  belongs  to  the  family  of  trapped 
acoustic  modes  propagating  in  a  wave-guide  between  the  wall  and  the  sonic  line  (Mack,  1984; 
Guschin  &  Fedorov,  1989).  The  existence  of  the  second  mode  was  established  by  the 
experiments  of  Kendall  (1975),  Demetriades  (1974),  Stetson  et  al.  (1983,  1992),  and  Kimmel  et 
al.  (1995).  Once  the  second  mode  sets  in,  it  becomes  the  dominant  instability  since  its  growth 
rate  tends  to  exceed  that  of  the  first  mode.  For  insulated  surfaces,  this  occurs  for  Mach  number 
larger  than  4.  For  cooled  surfaces,  the  second  mode  can  dominate  at  even  lower  Mach  numbers 
(Lysenko  &  Maslov,  1984). 

In  high-speed  flows,  the  second  mode  is  associated  with  disturbances  of  relatively  high  fi-equency 
corresponding  to  the  ultrasonic  band.  Malmuth  et  al.  (1998)  assumed  that  a  passive  ultrasonically 
absorptive  coating  (UAC)  of  fine  porosity  may  suppress  these  fluctuations  and,  at  the  same  time, 
may  not  trip  the  boundary  layer  due  to  roughness  effects;  i.e.,  the  passive  UAC  may  stabilize  the 
second  mode  by  a  disturbance  energy  extraction  mechanism.  This  hypothesis  was  examined  by 
the  inviscid  linear  stability  analysis.  Using  the  WKB  method,  Malmuth  et  al.  (1998)  formulated 
the  boundary  condition  on  an  ultrasonically  absorptive  wall  for  the  second  and  higher  modes,  and 
showed  that  the  absorption  does  cause  a  strong  stabilization  of  the  second  mode.  Later,  Fedorov 
&  Malmuth  (2001)  analyzed  the  absorption  effect  at  finite  Reynolds  numbers  using  viscous 
linear  stability  theoiy  and  found  that  viscosity  weakly  affects  the  stabilization  mechanism.  They 
also  considered  an  ultrasonically  absorptive  surface  of  a  particular  type;  namely,  the  wall  covered 
by  a  porous  coating  with  cylindrical  blind  micro-holes,  and  showed  that  a  relatively  thin  coating 
(of  the  thickness  ~0.1  of  the  boundary- layer  displacement  thickness)  can  dramatically  reduce  the 
second-mode  growth  rate. 

These  theoretical  findings  led  to  the  expectation  that  a  passive  porous  coating  may  be  exploited 
for  hypersonic  laminar  flow  control.  The  concept  was  verified  in  the  GALCIT  T-5  shock  tunnel 
of  Caltech  (Rasheed  et  al.,  2001)  by  testing  a  5”  half-angle  sharp  cone.  The  cone  has  a 
longitudinal  half  of  its  surface  solid  and  the  other  a  porous  sheet  perforated  with  equally  spaced 
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blind  cylindrical  holes  of  the  average  pore  diameter  60  pm,  depth  500  pm  and  average  spacing 
100  pm.  The  model  was  instrumented  by  thermocouples,  and  the  transition  onset  point  was 
determined  from  the  Stanton  number  distributions  measured  simultaneously  on  both  sides  of  the 
model  for  each  run.  The  experiments  were  performed  for  the  ranges  of  the  freestream  total 
enthalpy  4.18  <//(,  <13.34  MJ/kg  and  the  freestream  Mach  number  4.59  <M„  <6.4.  The 

porous  coating  delayed  transition  by  a  significant  amount.  For  the  majority  of  runs,  the  boundary 
layer  on  the  porous  surface  was  laminar  up  to  the  model  base,  while  transition  on  the  solid 
surface  was  observed  at  a  half  of  the  cone  length.  The  experiments  qualitatively  confirmed  the 
theoretical  prediction  of  Fedorov  &  Malmuth  (2001).  However,  quantitative  comparison  was  not 
feasible  because  the  cone  was  not  long  enough  to  measure  the  transition  locus  on  the  porous 
surface.  Since  the  boundary-layer  disturbances  were  not  measured,  these  experiments  did  not 
give  direct  evidence  of  the  second-mode  instability,  and  its  effect  on  the  transition  process  was 
not  clear.  This  motivated  us  to  perform  a  series  of  stability  experiments  on  sharp  cones  covered 
by  porous  coatings  of  various  microstructures.  Further  theoretical  research  was  also  needed  to 
address  randomly  porous  coatings  as  well  as  to  account  for  rarefied  gas  effects  on  the  UAC 
performance. 

The  first  series  of  stability  experiments  (Fedorov  et  al,  2003 a,b)  was  conducted  in  the  Mach  6 
wind  tunnel  of  the  Institute  Theoretical  of  Applied  Mechanics  (ITAM,  Novosibirsk,  Russia)  on  a 
7“  half-angle  sharp  cone  whose  longitudinal  half  surface  was  solid  and  other  half  surface  was 
covered  by  a  thin  porous  coating  of  random  structure,  namely,  a  fibrous  absorbent  material  (felt 
metal).  Hot-wire  measurements  of  i  natural!  disturbances  and  artificially  excited  wave  packets 
were  performed  on  both  solid  and  porous  surfaces.  Stability  analyses  for  2-D  and  3-D 
disturbances  showed  that  the  porous  coating  stabilizes  the  second  mode  and  marginally 
destabilizes  the  first  mode.  These  results  are  in  a  qualitative  agreement  with  the  experimental 
data  for  natural  disturbances.  The  theoretical  predictions  (including  nonparallel  effects)  are  in  a 
good  quantitative  agreement  with  the  stability  measurements  for  artificially  excited  wave  packets 
associated  with  the  second  mode  (Fedorov  et  al.,  2003a,b). 

In  this  report,  we  discuss  stability  experiments  and  calculations,  which  has  been  performed  for  a 
7°°half-angle  sharp  cone  with  the  porous  coating  of  regular  microstructure  similar  to  that  tested  in 
the  GALCIT  T-5  shock  tunnel  (Rasheed  et  al.,  2001).  In  Chapter  I  we:  formulate  stability 
problem  for  disturbances  in  the  boundary  layer  on  porous  wall;  investigate  propagation  of 
disturbances  within  pores  including  the  gas  rarefaction  effect  and  specify  boundary  conditions  for 
the  boundary-layer  disturbances  on  the  UAC  surface;  perform  parametric  studies  of  the  UAC 
stabilization  effect  on  the  second  mode.  In  Chapter  n  we  present  benchmark  measurements  of 
UAC  characteristics  and  their  comparisons  with  the  theoretical  model.  In  Chapter  HI  we:  discuss 
exj)erimental  setup,  a  sharp  cone  model  and  instrumentation;  describe  the  measurement 
techniques  and  data  processing;  present  results  of  stability  measurements  for  natural  and  artificial 
disturbances,  compare  experimental  data  with  predictions  of  linear  stability  theory.  Then  we 
conclude  the  report  and  discuss  future  efforts. 
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Chapter  I.  Theoretical  modeling  of  UAC  Performance 
1 .1 .  Linear  Stability  Analysis 


1.1.1.  Mean  flow 

For  the  laminar  boundary  layer  on  a  flat  plate  or  a  sharp  cone  at  zero  angle  of  attack,  equations 
are  written  in  the  self-similar  form  (Hayes  &  Probstein,  1959) 

(Cf)+ff''  =  0,  (1-1) 


+  /5'4 


U' 


*2 


h: 


/ 

/'/" 


=  0, 


(1.2) 


C  = 


PJi_ 

Pith 


H\ 


*2 


(7-m' 


^  =  |14 
P 


(1.3) 


Here  f{ri)  =  U*  /U],  g{ri)  =  H*  / ,  p,  /x,  Pr  and  7  are  the  streamwise  velocity,  total 
enthalpy,  density,  viscosity,  Prandtl  number  and  specific  heat  ratio,  respectively;  primes  denote 
differentiation  with  respect  to  17;  asterisk  stands  for  dimensional  quantities.  The  independent 
variables  are  given  by  the  Howarth-Dorodnitsyn  transformations  and  Mangier  transformation  of 
coordinates 


(1.4) 


where  x  is  coordinate  along  an  external  streamline;  y  is  coordinate  normal  to  the  body  surface; 

is  distance  from  the  symmetry  axis  to  the  wall  surface;  j  =  0  for  planar  flow  and  j  =  1  for 
axisymmetric  flow.  It  is  assumed  that  the  porous  coating  weakly  affects  the  mean  flow  weakly, 
and  the  standard  boundary  conditions  are  imposed  on  the  wall  surface 

T]  —  0:  /  =  /'  =  0 ,  5  =  (or  p'  =  0  for  adiabatic  wall) ,  (1.5) 

?7— >00:  f  =  1,  g  =  1.  (1-6) 


In  all  calculations  discussed  hereafter,  the  fluid  is  a  perfect  gas  of  7  =  1.4  and  Pr  =  0.708 .  The 
viscosity-temperature  dependency  is  approximated  by  Sutherlandis  law 


fi{T)  = 


(1  +  y3/2 

(T  +  S) 


(1.7) 


where  5  =  110/T/;  p,  =  p  /  p]  and  T  =  T*/T*  are  nondimensional  viscosity  and 
temperature. 
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1.1.2.  Linear  stability  problem 

Stability  analysis  includes  nonparallel  effects  associated  with  downstream  growth  of  the 
boundary-layer  thickness.  Effects  relevant  to  the  cone  surface  curvature  and  conical  divergence 
of  streamlines  are  neglected.  The  analysis  is  based  on  the  method  of  multiple  scales  applied  by 
Gaster  (1974)  to  incompressible  boundary  layers  and  extended  by  Padhye  &  Nayfeh  (1979)  to 
compressible  flows.  This  method  (with  insignificant  variations)  has  been  used  to  investigate  the 
nonparallel  effect  on  stability  of  supersonic  boundary  layers  (Gaponov,  1980;  El-Hady,  1980; 
Tumin  &  Fedorov,  1982).  Its  extension  to  the  case  of  multiple  modes  was  discussed  by  Zhigulev 
&  Tumin  (1982)  and  Fedorov  &  Khokhlov  (2002). 

Following  Fedorov  &  Khokhlov  (2002)  we  outline  the  nonparallel  linear  stability  analysis  for  a 
monochromatic  wave.  The  coordinates  {x,y,z)  are  made  nondimensional  using  the  boundary- 

layer  scale  r  =  /U]  ,  where  the  distance  V  from  the  leading  edge  is  assumed  to  be  much 

larger  than  T,  and  the  ratio  e  —  T/L*  is  treated  as  a  small  parameter.  Time  t  and  pressure  P 
are  referenced  to  I*  /Ul  and  respectively;  other  flow  characteristics  are  made 

nondimensional  using  their  quantities  at  the  upper  boundary-layer  edge.  Introducing  the  slow 
variable,  =  x*  / E  =  ex,  we  specify  the  mean-flow  velocity  components  {U,  V)  and 

temperature  T  as 

U  =  U{x^,y),  V  =  eV^(x^,y),  T  =  T{x^,y) .  (1.8) 


A  monochromatic  disturbance  is  represented  by  the  vector  function 


.  du  „  69 

Z  =  (u,  — ,  V,  p,  0,  — ,  w, 
dy  dy 


dw^T 


(1.9) 


7>{x, y, z,  t)  =  F(a;, y) exp{il3z  -  iut) , 


where  u,  v,  w,  p  and  9  are  velocity  components,  pressure  and  temperature;  jS  =  /?*/*  and 
ijj  =  ij/l*  /Ul.  The  amplitude  vector-function  'P{x,y)  satisfies  a  system  of  partial  differential 
equations  that  resulted  from  Fourier  transform  of  the  linearized  Navier-Stokes  equations  with 
respect  to  time  and  z -coordinate.  These  equations  with  the  boundary  conditions  are  written  in 
the  matrix-operator  form 


y  =  0:  f;  =  f;  ==  f;  =  f;  = 

y^oo:  Fj,  F;,  F^,  F^  -^0. 


(i.id) 

(1.11) 

(1.12) 


Here  the  admittance  components  A^,  A^,  A^  and  the  thermal  admittance  A^  are  complex 

quantities,  which  couple  velocity  and  temperature  disturbances  with  the  pressure  disturbance  on 
the  porous  surface.  The  admittances  depend  on  properties  of  the  porous  material,  mean-flow 
characteristics  on  the  wall  surface,  and  disturbance  frequency.  These  dependencies  are  specified 
in  Section  1.2  for  the  UAC  of  regular  microstructure.  On  the  solid  wall 
A^  —  All  —  A^  —  A^  =  0,  which  corresponds  to  no-slip  conditions  and  zero  temperature 


14 


perturbation  on  the  surface  of  high  thermal  conductivity. 


We  consider  a  partial  solution  of  the  problem  (1.10)-(1.12),  which  is  a  discrete  mode  of  the 
complex  eigenvalue  a{x^,P,u}) .  The  amplitude  vector-function  is  expressed  in  the  form 


F  =  (Fj  4- eFj  +  ...)exp{ie  ^S),  S  =  J a(x^)dXj^ . 


(1.13) 


Substituting  (1.13)  into  (1.10)-(1.12)  and  grouping  terms  of  the  same  order  of  magnitude  with 
respect  to  e  we  obtain  a  sequence  of  problems  for  F^(x^,p),  j  =  0,  1,  ...  In  the  zero-order 

approximation,  the  eigenvalue  problem  is  written  as 


1^2/ 


0 


(1.14) 


y  =  0 :  =  4^;, ,  ^3  =  44^ ,  ^3  =  44, ,  4^  =  44, ,  (i.i 5) 

2/^00:  4i,  4,3,  4^,  47^0.  (1.16) 


This  problem  describes  stability  of  a  locally  parallel  mean  flow.  Here  the  matrix  has 

dimension  8x8;  its  non-zero  elements  are  given  in  Appendix.  Solution  of  the  problem  (1.14)- 
(1.16)  is  expressed  in  the  form 

*1)  =c(a:i)C(a:i,y,«).  (l-l^) 


where  C  is  eigenvector  normalized  by  the  condition  5  i  ®-»  pressure- 

disturbance  amplitude  on  the  wall  surface  is  c(a:J.  For  spatial  instability  of  two-dimensional 
boundary  layers,  the  angular  frequency  u>  and  the  transverse  wavenumber  component  /?  are  real, 
whereas  a  is  a  complex  eigenvalue.  If  Im  a  <  0 ,  then  the  disturbance  amplifies  downstream 
with  the  spatial  growth  rate  a  —  —a..  The  amplitude  coefficient  c(x^)  is  determined  from  the 

first-order  approximation,  which  leads  to  the  inhomogeneous  problem 


fA_H 


F  --i^^  +  HF 
'  da  dx,  '  “ 


y  0  •  4l  ’  ^03  A'^04  ’  ^05  ^9^04  ’  '^7  ’ 


(1.18) 

(1.19) 


y— >oo:  4,)  43> 


07 


(1.20) 


The  right-hand  side  of  (1.18)  is  associated  with  the  nonparallel  effect;  the  vector  G  =  H^F,,  is 

given  in  Appendix.  The  problem  (1.18)-(1.20)  has  a  nontrivial  solution  if  its  right-hand  side  is 
orthogonal  to  the  eigenvector  ^(x^,y,a)  of  the  adjoint  problem 


(1.22) 


t/  =  0  :  ^2  =  Ig  =  4  =  0  ,  ^4  +  4r4  +  A^3  +  A4  +  ^^^7  “ 

y  ^oo:  4,  4,  4,  4 -"0,  (1-23) 

where  the  overbar  denotes  a  complex  conjugate  value.  The  orthogonality  condition  leads  to  the 
ordinary  differential  equation  for  the  amplitude  function  c{x^) 


4^C 


dc 

dx^ 


’  da  dx^l 


(1.24) 


where  the  scalar  product  is  defined  as 

00  g 

(1-25) 

0  7=1 

Substituting  solution  of  (1.24)  into  (1.17)  and  (1.13),  we  express  the  amplitude  vector  in  the  form 
F  =  lCgC(x^,y)  +  C)(e)]exp(ie“^  J[a;(Xj)  +  £W{x^)  +  0{e^)]ix^y  (1.26) 


where  is  constant.  The  amplification  rate  of  any  physical  quantity  can  be  calculated  using 
(1 .26)  and  (1 .27).  For  example,  the  x  -component  of  mass  flow  disturbance  is 

Q{x^,y)  =  [c^q{x^,y)  +  0(sr)]exp|ie“^  J[a;(2:j)  +  £W{x^)  +  C>(e^)]da;j) ,  (1.28) 


q  =  K+{lMX-^jT)U]lT.  (1.29) 

The  disturbance  growth  rate  is  a  logarithmic  derivative  of  (1.28),  which  is  expressed  in  the  form 

{d\a\q{x^,y)\ 


a{x^,y)  =  -a.  +  e 


dx. 


-W. 


+  0i£^). 


(1.30) 


Due  to  the  nonparallel  effect,  the  growth  rate  depends  on  the  vertical  coordinate  y.  In 
accordance  with  the  experimental  measurements  of  Chapter  in,  cr{x^,y)  is  calculated  at  y  —  y^ 

relevant  to  the  maximum  of  mass-flow  disturbance  in  the  boundary  layer.  This  maximum  lies  in 
the  critical  layer,  where  the  mean-flow  velocity  is  close  to  the  phase  velocity. 

The  zero-order  and  first-order  approximations  of  a  were  validated  by  comparison  with  the 
experimental  data  of  Kendall  (1967)  and  the  calculations  of  Chang  et  a/.  (1991)  performed  in  the 
framework  of  linear  parabolized  stability  equations  (PSE).  Figure  I.l  shows  the  first-mode 
growth  rate  as  a  function  of  the  nondimensional  frequency,  F  =  wVJ  /  ,  for  the  boundary 
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layer  on  a  flat  plate  at  the  Mach  number  =  2.2  and  the  Reynolds  number 

R  =  I  v\  =  1000.  The  wave  angle  is  ^  =  arctan(/3/Q;^)  =  60°  that  approximately 

corresponds  to  most  unstable  waves.  The  results  are  in  good  agreement  with  the  linear  PSE 
calculations  that  confirms  the  method  discussed  above. 


Fx10® 


Fig.  LI.  The  first-mode  growth  rate  as  a  function  of  frequency;  boundary  layer  on  a  flat  plate  at 
M^=2.2,  i?=1000,  the  wave  angle  }/=  arctan(yS/cirJ  =  60° :  triangles,  the  experiment  of 

Kendall  (1967);  dashed  line,  parallel  theory  (pres.);  solid  circles,  parallel  theoiy,  Chang  et  al. 
(1991);  solid  line,  non-parallel  theory  (pres.);  solid  squares,  non-parallel  PSE,  Chang  et  al. 
(1991). 


1 .2.  Acoustic  properties  of  UAC 

In  this  section,  the  solution  describing  propagation  of  monochromatic  acoustic  waves  inside  a 
cylindrical  pore  filled  with  low-density  gas  is  obtained.  The  analysis  is  focused  on  the  most 
practical  case,  when  corpuscular  effects  are  appreciable  only  inside  thin  Knudsen  layers  forming 
on  solid  surfaces.  Beyond  Knudsen  layers  gas  is  treated  as  a  continuum  media.  Its  dynamics  is 
governed  by  the  linearized  Navier-Stokes  equations  with  the  boundary  conditions  of  the  first- 
order  approximation  with  respect  to  Knudsen  number.  The  boundary  conditions  include  a 
nonzero  relative  macroscopic  gas  velocity  and  temperature  jumps  on  the  pore  surface.  Acoustic 
properties  of  the  perforated  panel  tested  in  the  wind  tunnel  T-326  (see  Chapter  III)  are  discussed. 
It  is  shown  that  the  rarefied  gas  effects  can  lead  to  essential  changes  of  acoustic  absorption  in  the 
range  of  nondimensional  parameters  relevant  to  the  wind-tunnel  conditions. 
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1.2.1.  Background 


Solutions  of  Navier-Stokes  equations,  which  account  for  effects  of  energy  dissipation  in  acoustic 
waves  traveling  inside  a  circular  tube,  are  known  since  the  latter  half  of  XIX  century  (Kirchoff, 
1868;  Strutt,  1894).  These  solutions  are  widely  used  in  analyses  of  acoustical  properties  of  wave¬ 
guides  and  perforated  panels  as  well  as  in  developments  of  semi-empirical  models  predicting 
sound  absorption  by  rigid-frame  porous  media  (Zwikker  &  Kosten,  1949;  Biot,  1956a,b;  Delany 
&  Bazley,  1970;  Attenborough,  1987;  Johnson  et  al,  1987;  Allard  &  Champoux,  1992). 

The  flow  resistance  and  heat  conductivity  of  micro-porous  materials  strongly  depend  on  ambient 
pressure  (Deryagin,  1946;  Dulnev  et  al.,  1974).  This  dependency  is  associated  with  molecular 
nature  of  gas  inside  the  pores.  The  interaction  of  gas  molecules  with  solid  surfaces  and  molecular 
collisions  in  space  obey  different  laws.  This  leads  to  formation  of  Knudsen  layers  on  solid 
surfaces.  The  Knudsen  layer  thickness  is  of  the  order  of  molecular  mean  free  path  A .  Beyond  the 
Knudsen  layer,  gas  can  be  treated  as  a  continuum  media.  However,  molecular  processes  inside 
the  BCnudsen  layer  affect  the  boundary  conditions.  On  macroscopic  level,  the  rarefaction  effects 
lead  to  slipping  and  temperature  jumps  on  the  interface  boundaries. 

Kundt  &  Warburg  (1875)  were  the  first  who  revealed  experimentally  the  slipping  in  the  gas 
outflow  from  a  circular  tube.  They  showed  that  the  volume  rate  of  laminar  low-density  gas  flow 
inside  the  tube  exceeded  the  volume  rate  predicted  by  the  Poiseuille  law.  Maxwell  (1952)  gave 
qualitative  explanation  to  this  phenomenon  by  means  of  molecular-kinetic  consideration. 
Smoluchowski  (1899)  discovered  physical  reasons  causing  the  temperature  discontinuity  on  the 
surface  of  a  solid  body  submerged  into  non-uniformly  heated  gas  of  low  density.  Since  the  work 
of  Kundt  &  Warburg  (1875),  the  mentioned  above  effects  were  thoroughly  investigated  using 
experimental  and  theoretical  methods.  Detail  reviews  of  these  studies  are  given  in  (Goodman  & 
Wachman,  1976;  Patterson,  1956;  Shidlovski,  1965)  and  many  other  monographs. 

Since  the  molecular  mean  free  path  A  is  inversely  proportional  to  the  pressure,  the  molecular 
effects  become  appreciable  in  the  whole  gas  volume  at  sufficiently  low  pressures,  and  the 
described  above  model  is  invalid.  However,  in  this  case  a  contribution  of  the  rarefied  gas  into 
averaged  macroscopic  properties  of  a  porous  material  is,  as  a  rule,  so  small  that  its  theoretical 
modeling  is  not  interesting  for  practice. 

Slipping  and  temperature  discontinuities  on  the  interface  boundaries  are  observed  in  media  with 
non-uniform  distributions  of  velocity  and  temperature.  Such  non-uniformities  are  continuously 
generated  by  acoustic  waves  propagating  within  pores.  As  a  consequence,  acoustic  properties  of 
porous  media  should  be  exposed  to  the  molecular  effects  at  sufficiently  low  pressures.  For  the  T- 
326  wind  tunnel  conditions,  the  Knudsen  number  is  Kn  =  X/r^  ~  0.1  (r^  is  pore  radius  of  the 

perforated  sheet),  and  acoustic  absorption  of  the  porous  layer  may  be  quite  different  from  that 
predicted  in  the  framework  of  continuum  fluid.  Similar  situations  may  occur  in  many  practical 
cases  such  as  propagation  of  sound  in  micro  ducts.  However,  we  have  not  found  any  publications 
on  this  subject.  Hereafter  we  discuss  results  of  our  theoretical  and  experimental  studies  of  the 
rarefaction  effects  on  acoustic  properties  of  the  porous  material  comprising  isolated  cylindrical 
micro-pores. 

1.2.2.  Problem  formulation 

Consider  propagation  of  a  monochromatic  acoustic  wave  inside  an  infinitely  long  circular  tube  of 
radius  ,  which  represents  an  isolated  pore.  The  tube  is  filled  with  gas  of  low  density.  The 

molecular  effects  are  assumed  to  be  appreciable  in  the  Knudsen  layer  only.  Beyond  the  Knudsen 
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layer  and  far  from  the  tube  ends,  the  disturbance  generated  by  infinitesimal  oscillating  pressure 
gradient  Re|(dpy/d2:)exp(ia;t)|  is  governed  by  the  linearized  momentum  and  energy  equations. 

For  the  amplitudes  of  velocity  u  and  temperature  f ,  these  equations  can  be  written  in  the  form 


Au  4-  A^u  =  —— 
fj,  dz 


AT  +  PrA'f  =  -  — p  , 

K 


(1.31) 


where  the  z  -axis  is  directed  along  the  pore  axis,  the  origin  of  the  reference  system  is  placed  in 
the  area  center  of  the  tube  cross-section;  i  =  u  =  27r/  is  circular  frequency;  Pr  is  Prandtl 
number  (Pr  =  k);  is  specific  heat  capacity  under  constant  pressure;  p  and  k  are 

dynamic  viscosity  and  heat  conductivity  of  gas  inside  the  tube;  A  =  ^—uop^r^j^ ;  is  gas 


1  d 

density  in  equilibrium;  A  —  — — - 

r  or 


•A' 

dr  j 


is  Laplacian,  where  r  is  radius  dimensionalized  using 


the  pore  radius  r^. 


In  the  first-order  approximation  with  respect  to  small  (but  finite)  Knudsen  number,  the  boundary 
conditions  are 


“s  =  -AKn 


(1.32a) 


2;  = 


KnidT' 

^  Priori’ 


(1.32b) 


where  and  are  gas  velocity  and  temperature  on  the  solid  surface;  the  dimensionless  factors 
and  depend,  mainly,  on  laws  of  interaction  between  gas  molecules  and  the  solid  surface; 
the  subscript  i  wi  denotes  quantities  on  the  pore  wall. 

In  the  framework  of  elementary  kinetic  theory  (see,  for  example,  Goodman  &  Wachman,  1976), 
the  dimensionless  coefficients  A^ ,  A^  and  the  BCjiudsen  number  Kn  can  be  expressed  in  terms 

of  the  molecular  tangential  impulse  and  energy  accommodation  coefficients  and  other 
macroscopic  parameters  as 

4,  =  -  0.5,  A^  =  27  (a-^  -  0.5)/(7  + 1) , 


Kn  =  4p/(p„crJ, 


where  7  is  specific  heat  ratio;  c  =  is  mean  module  of  molecular  velocity  vector; 

is  gas  constant  per  unit  mass;  is  gas  temperature  in  equilibrium.  In  addition,  we  assume  that 
the  solid  surface  temperature  in  (1.32b)  is  constant.  The  boundary  conditions  for  more 
general  cases,  in  particular,  for  rapidly  varying  wall  temperature  T^,  are  given  by  Patterson 
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(1956)  and  Shidlovski  (1965). 

The  discussed  above  formulation  is  based  on  the  assumption  that  the  Knudsen  layer  thickness  is 
much  less  than  the  pore  radius  (A  <  r^).  This  allows  us  to  use  the  continuum  media  model 

beyond  the  Knudsen  layers.  However,  besides  the  Knudsen  number  the  problem  depends  on  the 
ratio  of  the  mean  free  path  to  the  Stokes  layer  thickness  6^  ~  ■  As  a  result,  the  solution  is 

a  function  of  both  Kn  and  the  product  Kn  | A| ,  which  should  be  small  also.  The  condition 
Kn|A|  ■C  1  becomes  important  for  disturbances  of  relatively  high  frequency,  when  the  Stokes 
layer  thickness  is  less  (or  much  less)  than  the  pore  radius;  i.e.,  |A|  >  1 . 


1.2.3.  Solution  of  the  problem 

Using  the  substitutions,  which  are  valid  for  the  time  factor  exp{iu>t) , 

u[r,z)  =  -  F{r,B^,A)}, 

lup^  dz  - 


(1.33) 

(1.34) 


where  =  AJAn,  =  (A^/PrjKn,  we  obtain  that  the  function  F(r,BX)  is  a  solution  of 
the  Bessel  equation 


[r  dr  dr 


F{r,BX)  =  0, 


(1.35) 


with  the  boundary  condition 

dr 


(1.36) 


lr=l 


Here  B  =  B^  and  ^  =  A  in  the  expression  (1.33)  for  the  velocity  amplitude;  B  =  B^  and 
C  =  A^^  in  the  expression  (1.34)  for  the  temperature  amplitude. 

The  solution  of  (1 .35),  (1.36)  is  expressed  it  terms  of  Bessel  function  of  first  kind  of  zero  order 

F{r,B,0  =  CJ^{Cr),  (1-37) 

where  C  is  integration  constant,  which  is  determined  from  the  boundary  condition  (1.36)  as 

1 


C  = 


T(C)[1-0.55C%)] 


C^n(C) 


(1.38) 


1 

Averaging  (1.33)-(1.34)  over  the  pore  cross-section  as  obtain 
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(1.39) 

(1.40) 

Finally,  using  the  general  relations 

wp{u){u{z))  = 

(1.41) 

^  1  (^)  iL 

(1.42) 

where  is  unperturbed  ambient  pressure,  we  can  express  the  complex  dynamic  density  p(u) 

and  the  dynamic  compressibility  C{uj)  as 

p(u)  =  pJ{l-{F{B„A))}. 

(1.43) 

C(a,)  =  -i-{l  +  (7-l)(F(B,,7P?A))}  , 

7Po 

(1.44) 

1-0.5SC%)‘ 

(1.45) 

In  the  case  of  Kn  =  0 ,  the  formulas  (1.43)-(1.45)  coincide  with  (2.3),  (2.4)  discussed  in  Chapter 

n. 


1.2.4.  Numerical  examples 

The  absorptive  coating  tested  in  the  benchmark  experiments  is  a  panel  with  equally  spaced 
cylindrical  pores  of  circular  cross-section.  The  mean  pore  diameter  is  60  |xm,  the  pore  depth  is 
500  pm,  and  the  mean  spacing  between  pores  is  100  pm.  Schematic  picture  of  this  perforated 
panel  is  shown  in  Fig.  1.2. 

The  theoretical  expressions  (2.43)-(2.45)  for  dynamic  density  and  dynamic  compressibility  are 
used  for  evaluation  of  the  ambient  pressure  effect  on  acoustic  properties  of  the  cylindrical  pore. 
The  Bessel  functions  of  complex  argument  are  calculated  using  the  approximations  of  the  Kelvin 
functions  presented  in  (Abramowitz,  &  Stegun,  1964). 

Calculations  were  conducted  for  pressures  =1.013  and  10.13  kPa  at  the  gas  temperature 
Tp  =  288.16  K.  The  disturbance  frequency  varied  from  0.8  kHz  to  10  kHz.  The  molecular 
tangential  impulse  (a^)  and  energy  (a^)  accommodation  coefficients  are  0 

the  pressure  1.013  kPa,  the  Knudsen  number  is  Kn  =  A/rQ=4.2- 10“\  For  the  majority  of 
calculations,  |A|  <  1  corresponds  to  low-frequency  disturbances. 
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Fig.  1.2.  Schematic  of  the  wall  covered  by  the  perforated  sheet. 


Results  of  calculations  are  shown  in  Figs.  I.3-L5.  The  real  part  of  dynamic  density,  Re(p) ,  is 
practically  constant  in  the  considered  frequency  band  and  fixed  pressure  (see  Fig.  1.3).  The 
rarefaction  effect  leads  to  essential  decrease  of  Re(p)  at  =  1.013  kPa.  For  the  ambient 

pressure  10.13  kPa,  this  effect  is  negligible.  Figure  1.4  shows  that  the  imaginary  part  of  dynamic 
density,  Im(p),  is  more  sensitive  to  pressure  variations  than  its  real  part.  This  sensitivity 
increases  as  the  disturbance  frequency  decreases.  For  sound  of  low  frequency  (/  w  1  kHz),  the 
rarefaction  effect  leads  to  approximately  two  times  increase  of  Im(p)  at  =  1.013  kPa. 

The  real  part  of  dynamic  compressibility  is  practically  constant,  Re|C'(a;)7Pg}  «  1.4,  in  the 
considered  ranges  of  pressure  and  frequency.  Figure  1.5  shows  that  the  imaginary  part, 
Im|C(a;)7Pg} ,  is  approximately  a  linear  function  of  frequency.  The  slope  of  this  function 
essentially  depends  on  the  ambient  pressure  and  Knudsen  number.  However,  the  module  of 
Im|(7(a;)7Pp}  is  much  smaller  than  Re|C'(w)7Pg}  and  can  be  neglected  in  predictions  of 

acoustic  properties  of  the  perforated  panel. 
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Fig.  1.3.  Theoretical  dependence  of  the  real  part  of  dimensionless  dynamic  density  for  the 
cylindrical  pore  of  radius  30p,m  versus  pressure  and  frequency  (red  solid  and  red  dashed  lines 
coincide). 


Fig.  1.4.  Theoretical  dependence  of  the  imaginary  part  of  the  dimensionless  dynamic  density  for 
the  cylinder  pore  of  radius  30pm  versus  pressure  and  frequency. 
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Fig.  1.5.  Theoretical  dependence  of  the  imaginary  part  of  the  dimensionless  dynamic 
compressibility  for  the  cylinder  pore  of  radius  30pm  versus  pressure  and  frequency. 


1.2.5.  Admittance  of  porous  layer 

As  shown  by  Fedorov  &  Malmuth  (2001)  and  Fedorov  et  al.  (2003 a, b),  the  admittances 
and  A^  produce  small  effects  on  the  boundary-layer  instability,  compared  to  the  effect  of  A^ . 
This  allows  us  to  approximate  the  wall  boundary  conditions  (1 .1 1)  as 

f;(o)  =  o,  f;(o)  =  o,  f;(o)  =  o,  (i.46) 

f;(o)  =  4f;(o). 

The  admittance  component  A^  is  expressed  in  the  form  (Fedorov  &  Malmuth,  2001) 

A^  = —ta,nh.{mh) ,  (1-47) 


where  ^  is  porosity,  h  —  h*/l*  is  the  porous-layer  thickness,  and  m  are  the  characteristic 
impedance  and  propagation  constant  of  a  porous  medium,  respectively.  The  two  latter  are 
functions  of  the  complex  dynamic  density  p  and  the  complex  dynamic  compressibility  C  .  Using 
the  scaling  of  stability  problem  formulated  in  Section  1.1,  we  obtain 


^/?7g 


(1.48) 
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where  the  dimensionless  complex  dynamic  density  p 
expressed  in  the  form 

and  dynamic  compressibility  C  are 

(1.49) 

a  =  l  +  (7-l)F(5^,?^^), 

(1.50) 

(1.51) 

F(n  -Jf^\  Q(cVPr) 

i-0.5S,kVP?)'(3{sVpF)’ 

(1.52) 

Note  that  in  the  boundaiy-layer  stability  problem  the  time  factor  is  exp{-iu>t).  For  this  reason, 
the  argument  ,  which  is  used  in  stability  analyses,  is  different  from  the  argument 

A  —  ^-icjp^r^  in  the  relations  (1.43)-(1.45)  derived  for  disturbances  -exp(ia;t).  The 

formulas  (1.49)-(1.52)  coincide  with  those  of  Fedorov  &  Malmuth  (2001)  used  for  the  modeling 
of  at  Kn  =  0 . 


1.3.  Parametric  studies  of  the  UAC  effect  on  the  second  mode 

In  this  section,  we  discuss  parametric  studies  of  the  UAC  effect  on  the  second-mode  instability 
for  the  T-326  wind  tunnel  conditions.  These  studies  are  performed  in  the  framework  of  local- 
parallel  linear  stability  theory;  i.e.,  the  growth  rate  (1.30)  is  calculated  at  £  =  0 .  The  nonparallel 
effects  will  be  included  in  Chapter  HI  for  comparison  of  stability  calculations  with  the  wind- 
tunnel  data. 

1.3.1.  Basic  parameters  and  mean  flow 

Calculations  are  performed  for  a  7°  half-angle  sharp  cone  at  the  freestream  Mach  number 
w  6 .  Flow  parameters  at  the  upper  boundary  layer  edge  are: 

1.  Mach  number  =  5.3 

2.  Temperature  T*  =  56.4  K 

3.  Stagnation  temperature  calculated  at  the  specific  heat  ratio  7  =  1.4 

r;  =  (1  +  (7  -  i)M,'  /2)r;  =  373.255  k 
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4.  Unit  Reynolds  number  Re^  =  U^/  is  varied  in  the  range  Re^  =  (8.5  —  47.8)  x  10®  m  ’ 
The  cone  length  is  I-  =  0.5  m.  Its  surface  is  thermally  insulated:  T*  =  T*^ . 

The  porous  coating  characteristics  correspond  to  the  experiments  in  T-326  wind  tunnel: 

1 .  Average  pore  radius  r*  ==  25  pm 

2.  Coating  thickness  h*  is  variable 

3.  Average  pore  spacing  s*  =  100  pm 

4.  Porosity  (j)  ^  0.2 


Fig.  1.6.  Boundary-layer  displacement  thickness  at  various  unit  Reynolds  numbers  Re^; 
adiabatic  cone,  local  Mach  number  =5.3,  local  temperature  at  the  upper  boundary-layer 
edge  T,  =  56.4  K. 


The  eigenvalue  problem  (1.14)-(1.16)  is  solved  at  Prandtl  number  Pr  =  0.708  and  specific  heat 
ratio  7  =  1.4.  The  dynamic  viscosity  is  calculated  using  the  Sutherland  formula;  the  second 
viscosity  is  p*  =  0.8p* . 

Figure  1.6  shows  the  boundary-layer  displacement  thickness  6' {x’) .  The  nondimensional  mean- 
flow  profiles  U{y)  and  T{y) ,  y  =  y'  /  6* ,  are  shown  in  Figs.  1.7  and  1.8,  respectively.  The 

boundary-layer  thickness,  6  1.3S\  is  marked  by  arrow.  The  adiabatic  wall  temperature  is 

T  =5.707T. 

ad  ® 
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Fig.  1.7.  Mean-flow  velocity  profile;  adiabatic  cone,  local  Mach  number  =  5.3,  temperature 
at  the  upper  boundary-layer  edge  =  56.4  K. 


Fig.  1.8.  Mean  temperature  profile;  adiabatic  cone,  local  Mach  number  =5.3,  temperature  at 
the  upper  boimdary-layer  edge  =  56.4  K. 


1.3.2.  Stability  characteristics  for  adiabatic  wall 

We  consider  characteristics  of  the  second  mode.  Since  most  unstable  second-mode  waves  are 
two-dimensional,  all  calculations  are  performed  for  /5  =  0 .  Results  are  presented  for  the  unit 

Reynolds  numbers  Re^  =  plUH  i\  =  (10,  20,  30,  40)xl0®m'’  relevant  to  the  experimental 
conditions. 

Disturbances  of  maximum  growth  rate 

Figures  1.9  and  1. 10  show  the  disturbance  wavelength  A  =  27r/Re(Q;)  and  the  nondimensional 
frequency  parameter  F  —  /  Ul  for  the  second-mode  waves  of  maximum  growth  rates.  These 

characteristics  are  weakly  sensitive  to  the  porous  coating;  i.e.,  the  disturbance  wavelength  and 
frequency  on  the  porous  and  solid  walls  are  almost  equal.  The  ratio  X*  /S  k2  is  typical  for  the 
second-mode  instability. 

Figures  1.1 1-1.14  show  the  maximum  growth  rates  on  solid  and  porous  surfaces.  The  stabilization 
effect  increases  with  the  longitudinal  coordinate  X.  Maximum  stabilization  is  observed  at 
Re^  =  40  X 10”®  m'*  (see  Fig.  1.14)  that  corresponds  to  maximum  of  the  ratio  r*  /  6* . 

Figures  1. 15-1. 16  show  disturbance  characteristics  for  the  porous-coating  thickness 
h*  =2r,  4r*,  6r’,  8r*,  lOr*,  oo  at  the  unit  Reynolds  number  Re^  =  10  x  10®m‘’.  For 

h*  >  6r*,  the  reflection  from  the  pore  bottom  is  small;  i.e.,  the  porous  coating  can  be  treated  as 
infinitely  thick.  Similar  trends  are  observed  in  the  cases  of  Re^  =  (20  —  40)  x  10®  m'*. 


X,  mm 


Fig.  1.9.  Wavelength  of  the  most  unstable  second  mode  at  various  unit  Reynolds  numbers  Re^ ; 

T„  =  T^,M,  =  5.3,T,^56A  K. 
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X,  mm 


Fig.  1.10.  Frequency  parameter  of  the  most  unstable  second  mode  at  various  unit  Reynolds 
numbers  Re^;  =  T^,  =  5.3,  =  56.4  K. 
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Fig.  1.11.  Spatial  growth  rate  of  most  unstable  disturbance  on  solid  and  porous  surfaces; 
Re^  =  lO^rn-*,  r„  =  r,,  =  5.3,  =  56.4  K. 
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Fig.  1.12.  Spatial  growth  rate  of  the  most  unstable  disturbance  on  solid  and  porous  surfaces; 
Re^  =  2 X 10’’ m\  T^=T^,M^=  5.3,  T,  =  56.4  K. 


Fig.  1.13.  Spatial  growth  rate  of  the  most  unstable  disturbance  on  solid  and  porous  surfaces; 
Re ,  =  3  X  10^m-\  r„  =  T, ,  M,  =  5.3,  T,  =  56.4  K. 


X,  mm 


Fig.  1.14.  Spatial  growth  rate  of  the  most  unstable  disturbance  on  solid  and  porous  surfaces 
Re^  =  4xl0^m-',  r„  =  T,,  =  5.3,  =  56.4  K. 


Re,=10,E+6  1/m,  r=25  mkm,  s=100  mkm 


X,  mm 


Fig.  1.15.  Maximal  growth  rate  for  UAC  of  various  thickness;  Re^  =  lO’^m'*,  T„  =  , 

=  5.3,  =  56.4  K. 
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X,  mm 


Fig.  1.16.  Wavelength  of  most  unstable  waves  for  UAC  of  various  thickness;  Re^  =  10^  m  ', 
7;  =  7:.,M,  =  5.3,r,=56.4K. 

Disturbances  of  various  frequencies 

Figures  1.17  and  1.18  show  the  growth  rate  as  a  function  of  X  for  various  frequency  parameters 
F  =  u/v]  /  .  The  disturbance  of  fixed  frequency  is  stabilized  by  the  porous  coating  at  large 

X ,  while  it  is  destabilized  at  small  X .  The  net  effect  is  small  or  even  negative  for  relatively 
small  F.  This  behavior  is  typical  for  the  boundary  layer  on  adiabatic  wall.  Note  that 
aerodynamic  surfaces  of  actual  hypersonic  vehicles  correspond  to  small  wall-temperature  ratios: 
<0.3 .  In  this  case  the  UAC  performance  is  much  higher  as  will  be  shown  in  the  next 

section. 

1.3.3.  Stability  characteristics  for  cold  wall 
The  surface  temperature  =  32^  «  169K 

Figures  1.19  and  1.21  show  the  maximum  growth  rate  as  a  function  of  the  longitudinal  coordinate 
X  for  the  solid  wall  (black  line)  and  the  porous  wall  (red  line);  unit  Reynolds  number  is 
RBj  =  10^  and  2  x  lO’^  m'',  respectively.  Figures  1.20  and  1.22  show  the  growth  rate  as  a  function 

of  X  for  various  frequency  parameters  F  and  the  same  mean-flow  parameters.  In  this  case,  the 
net  stabilization  effect  is  positive. 

The  surface  temperature  =  27^  «  113  K 

Figures  1.23  and  1.25  show  the  maximum  growth  rate  as  a  function  of  the  longitudinal  coordinate 
X  for  the  solid  wall  (black  line)  and  the  porous  wall  (red  line);  unit  Reynolds  number 
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Re  =  10^  and  2  x  lO’'  m'',  respectively.  Figures  1.24  and  1.26  show  the  growth  rate  as  a  function 
of  X  for  various  frequency  parameters  F  .  In  this  case,  the  net  stabilization  effect  is  strong. 

These  calculations  indicate  that  the  UAC  performance  dramatically  increases  as  the  wall 
temperature  decreases.  This  trend  is  favorable  for  practice  because  surfaces  of  actual  hypersonic 
vehicle  have  wall  temperatures  essentially  lower  than  the  adiabatic  temperature. 


Re^=10.E+6  1/m,  M=5.3,  Tw=Tad,  r=25  mkm,  s=100  mkm,  Au=Ae=1,  Kn=0.79 


Fig.  1.17.  Second-mode  growth  rate  distributions  at  various  frequencies;  Rej=10^m'’, 


Re,=20.E+6  1/m,  M=5.3,  Tw=Tad,  r=25  mkm,  s=100  mkm,  Au=Ae=1,  Kn=0.395 


7  1 

Fig.  1.18.  Second-mode  growth  rate  distributions  at  various  frequencies;  Re^  =  2  x  10  m' , 
M,  =  5.3,  =  56.4  K. 
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Fig.  1.23.  Maximum  growth  rate  of  the  second  mode;  =  2T^ ,  Re^  =  10  m' , 
T,  =  56.4  K. 


Re,=10.E+6 1/m,  M=5.3,  Tw=2.,  r=25  mkm,  s=100  mkm,  Au=Ae=1,  Kn=0.188 
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Fig.  1.24.  Second-mode  growth  rates  for  various  frequencies;  =  2T^ ,  Re 
M,  =  5.3,  T  =  56.4  K. 
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1.3.4.  Rarefaction  effect 


Because  of  relatively  small  mean-flow  densities  the  rarefaction  effects  inside  pores  may  be 
important.  Hereafter,  the  Knudsen  number  is  determined  as 


Kn  = 


r’P: 


(1.53) 


where  =  287  J/(kg  K)  is  gas  constant.  The  Knudsen  number  can  be  also  expressed  in  terms  of 
basic  nondimensional  parameters  as 

=  (1.54) 

rR 

Figure  1.27  shows  the  disturbance  growth  rate  as  a  function  of  the  coating  thickness  h  with  and 
without  the  Knudsen  layer.  Calculations  were  performed  for  the  disturbance  of  fixed  frequency 
F  =  1.43  X 10”^  at  the  Reynolds  number  R  =  13262 .  The  rarefaction  effect  makes  pores  to  be 
more  transparent  with  respect  to  disturbances.  This  leads  to  stronger  stabilization  of  the  second 
mode. 

The  maximum  growth  rates  for  the  adiabatic  wall  covered  by  the  UAC  of  r*  =  25  pm, 
s*  —  100pm  and  various  h*  are  shown  in  Fig.  1.28  for  Re^  ==  10^ m''.  The  disturbance  reflection 
from  the  pore  bottom  becomes  negligible  for  h*  >  lOr* .  Similar  trends  were  obtained  for  the 
unit  Reynolds  numbers  Re^  =  (2  —  3)  x  10^  m'*. 


Data  shown  in  Figs.  I.29-I.31  were  calculated  for  infinitely  deep  holes  in  adiabatic  wall.  It  is  seen 
that  the  Knudsen  layer  on  pore  walls  leads  to  additional  stabilization  of  the  second  mode. 


T-326,  R=1.36258E+4,  Tw=4.84,  M=5.3,  F=1 .43047E-4,  r=0.02,  n=0.2 


Fig.  1.27.  Second-mode  growth  rate  as  a  function  of  the  UAC  thickness  with  and  without 
Knudsen  layer;  r„  =  4.84r, ,  R  =  13262 ,  F  =  1.43  x  10-“ ,  =  5.3,  =  56.4K. 
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T-326,  Re^=10.E+6  1/m,  r=25  mkm,  s=100  mkm,  Kn=0.79,  Au=Ae=1. 


1.4.  Summary 

Linear  stability  problem  for  three-dimensional  disturbances  in  two-dimensional  boundary  layers 
on  a  porous  wall  was  formulated  in  the  framework  of  the  asymptotic  method  of  multiple  scales. 
This  problem  differs  from  the  typical  stability  problem  because  of  new  boundary  conditions 
associated  with  absorption  of  disturbance  energy  by  a  porous  coating.  These  boundary  conditions 
are  formulated  in  terms  of  the  porous-wall  acoustic  admittances,  which  depend  on  the  UAC 
microstructure. 

Acoustic  properties  of  a  cylindrical  pore  were  analyzed  in  the  framework  of  linear  acoustic 
theory.  Dependencies  of  the  dynamic  density  and  dynamic  compressibility  on  gas  pressure  and 
disturbance  frequency  were  derived  including  the  gas  rarefaction  effects.  The  analysis  is  focused 
on  the  most  practical  case  when  the  molecular  mean  free  path  is  smaller  (or  much  smaller)  than 
the  pore  radius.  Calculations  indicate  that  acoustic  properties  of  the  perforated  panel  essentially 
depend  on  Knudsen  number  in  the  considered  ranges  of  pressure  and  frequency.  The  imaginary 
parts  of  dynamic  density  and  dynamic  compressibility  are  more  sensitive  to  the  pressure  variation 
than  their  real  parts. 

Parametric  calculations  of  the  UAC  effect  on  the  second-mode  stability  are  performed  for  the  T- 
326  wind  tunnel  conditions  including  the  gas  rarefaction  effect.  In  the  case  of  adiabatic  wall,  the 
net  stabilization  effect  is  small  or  even  negative  for  disturbances  of  relatively  small  frequencies. 
The  UAC  stabilization  dramatically  increases  as  the  wall  temperature  ratio  decreases.  This  trend 
is  preferable  for  practice  because  surfaces  of  actual  hypersonic  vehicle  are  maintained  at  low 
wall-temperature  ratios. 

The  rarefaction  effect  is  appreciable  for  the  T-326  conditions.  The  presence  of  the  Knudsen  layer 
on  the  pore  wall  leads  to  deeper  penetration  of  disturbances  into  pores  and  provides  stronger 
stabilization  of  the  second  mode. 

The  discussed  above  parametric  studies  were  used  for  design  the  porous  coating  tested  on  a  sharp 
cone  in  the  T-326  wind  tunnel  (see  Chapter  HI). 
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Chapter  II.  Benchmark  Measurements  of  UAC  Characteristics 

Design  of  ultrasonic  absorptive  coatings  (UAC)  for  hypersonic  boundary-layer  stabilization 
requires  sound-absorbing  characteristics  of  porous  materials.  It  is  difficult  to  perform  direct 
measurements  of  these  characteristics  in  aerodynamic  experiments  due  to  low  gas  pressure  (» 
0.01  bar)  and  high  frequencies  (/■»  10®  Hz)  of  disturbances  observed  in  hypersonic  boundary 
layers.  Under  these  conditions  acoustic  disturbances  quickly  attenuate  via  viscous  dissipation  of 
acoustic  energy.  The  attenuation  factor  is  proportional  to  f  and  inversely  proportional  to  the  gas 
density.  For  the  frequency  of  200  kHz  and  density  corresponding  to  the  pressure  of  20-30  ton- 
typical  for  free  stream  in  a  hypersonic  wind  tunnel,  sound  intensity  drops  e-times  at  the  distance 
of  few  millimeters. 

To  overcome  this  difficulty  it  is  suggested: 

•  Develop  a  theoretical  model  of  acoustic  absorption  by  a  porous  coating 

•  Perform  benchmark  measurements  at  low  frequencies  and  verify  the  theoretical  model 

•  Calculate  the  UAC  characteristics  for  required  ranges  of  frequencies  and  densities 
This  combined  experimental  and  theoretical  approach  is  accomplished  in  the  present  work. 

2.1.  Modeling  of  acoustic  propagation  in  porous  materials 

The  micrograph  of  the  porous  material  is  presented  in  Fig.  n.l  (the  large  scale  division  is  0.1 
mm).  Microscopic  measurements  showed  that  holes  have  conical  shape:  the  hole  diameter  is  =55 
micron  from  one  side  and  ^0  micron  from  the  other  side.  The  perforated  sheet  has  the  thickness 
^=0.5  mm  and  the  porosity  0  =0.3. 


Fig.  n.l  Micrograph  of  the  perforated  sheet  (top  view). 


The  porous  material  is  modeled  as  a  set  of  isolated  cylindrical  pores.  Stinson  &  Champoux 
(1992)  performed  theoretical  studies  of  acoustic  disturbances  propagating  in  cylindrical  pores  of 
arbitrary  cross-section.  The  disturbance  propagation  depends  on  viscous  effects  in  the  tube  and 
on  heat  conductivity.  These  effects  can  be  modeled  separately  (Stinson  &  Champoux,  1992; 
Zwikker  K.,  Kosten,  1952). 

Propagation  of  acoustic  disturbances  in  the  cylindrical  pore  is  characterized  by  the  wave 
resistance  W  and  the  propagation  constant  F.  These  values  can  be  expressed  in  terms  of  the 
complex  dynamic  density  p(G^f  and  the  complex  compressibility  C(ry)  as 


42 


W  =  [pirn)/ 

(2.1) 

V  =  iu{p(.u}')C  , 

(2.2) 

where  C  (w)  =  (l  /  7P0 )  {7  -  (7  - 1)  /  p  (Pr  a;)]| ,  which  is  valid  for  any  cylindrical  pore  of 

constant  cross-section;  po  is  air  density;  y  is  specific  heat  ratio;  Pr  is  Prandtl  number;  mis  angular 
frequency;  Po  n  air  static  pressure. 

For  a  cylindrical  pore  of  circular  cross-section,  the  wave  resistance  W  and  the  propagation 
constant  F  are  shown  in  Figs,  n.2  and  n.3  as  functions  of  the  nondimensional  parameter 
( 

T  P  U)  •  • 

H  —  -^  —  r  ,  where  is  pore  radius;  6^  is  acoustic  boundary-layer  thickness;  7]  is  an 

“I  ^  j 

dynamic  viscosity.  These  dependencies  were  calculated  using  the  relations  given  by  Benade 
(1968):  the  complex  density  is  expressed  as 

p  (w)  =  pjl  -  2T  ,  (2.3) 

where  T{x)  =  ix)/  and  J^-  Bessel  functions  of  zero  and  first  order;  the  complex 

compressibility  is 

C(^y)  =  (l//>oc'){l+(r-l)[2P(^/^^)/^/^^f]}  •  (2.4) 
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Fig.  n.2  The  wave  resistance  of  cylindrical  pore  of  circular  cross-section. 


Fig.  n.3  The  propagation  constant  of  cylindrical  pore  of  circular  cross  section. 


Consider  a  porous  coating  of  the  porosity  .  The  ratio  of  the  pore  length  to  the  coating  thickness 
is  q;  the  pore  tortuosity  is  defined  as  q^.  The  complex  density  of  the  porous  material  is  introduced 
by  the  relation 


wp^(u;)v^=-f,  (2.5) 

az 

where  is  the  mean  velocity  of  the  flow  through  the  total  cross  section  of  the  coating; 

—  {c^  /  .  The  coating  compressibility  is  determined  as  =  Pa)lP^ 

where  6^  is  acoustic  pressure  pulsations,  averaged  over  the  total  cross  section  of  the  coating; 
Cj  (fii)  =  ^  (fi?) .  Then  the  wave  resistance  and  the  propagation  constant  are  expressed  in  the 
form  similar  to  the  case  of  a  single  cylindrical  pore 

W  ,  (2.6) 

r  =  ia{pi  {co)C^  {(o))"^ .  (2.7) 

For  porous  materials  of  complex  microstructure,  details  of  the  pore  shape  are  not  available. 
Acoustic  characteristics  of  these  materials  are  determined  with  the  assumption  that  actual  pore 
structures  may  be  simulated  using  cylindrical  pores  of  effective  size.  For  this  purpose,  the 
following  non-acoustic  characteristics  of  the  porous  material  are  used:  porosity,  tortuosity,  and 
air  resistance  to  the  constant  stream.  These  parameters  can  be  determined  regardless  of  acoustic 
measurements.  In  addition,  so-called  shape  factor  is  introduced  as  a  scale  parameter  multiplier  to 
describe  the  material  via  cylindrical  pores.  In  this  framework,  asymptotic  approximations  for 
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various  frequency  bands  (namely,  low-frequency  band  for  small  H  and  high-frequency  band  for 
large  H)  are  normally  used. 

For  low  frequencies,  the  complex  density  and  compressibility  are  approximated  as 

pioj)  =  a  /  w  +  (2.8) 


C  (a;)  =  (l  /  7P0 )  {7  -  (7  - 1) \wp^  Pr /  cr]} , 


(2.9) 


where  \  is  the  structural  coefficient  independent  on  frequency  (it  is  determined  experimentally). 
According  to  Rzhevkin  (1960),  this  coefficient  is  called  as  i  added  mass  coefficient! . 


For  high  frequencies,  viscous  effects  are  appreciable  in  a  thin  layer  on  the  pore  wall  (acoustic 
boundary  layer).  In  this  case,  the  complex  density  and  compressibility  are 


piu})^  p^ 


l  +  (2/r„)(77/*a;po)'^' 


(2.10) 


C7(a;)  « 


1 


7P0 


l  +  (7-l)  — 

V 

0 


V 

1/2' 

(2.11) 


Stinson  &  Champoux  (1992)  performed  calculations  for  various  pore  cross  sections  and  showed 
that  the  frequency-independent  structural  factor  can  be  used  as  a  multiplier  to  the  scale  parameter 
to  find  the  dynamic  parameters  of  a  porous  material.  The  complex  density  is  expressed  in  the 
form 


pCu)  =  /Og  +  {a/iu})F{H). 
The  function  F  is  modified  as 


(2.12) 


1 


(2.13) 


The  scale  parameter  is  defined  as  H  =  ,  where  is  the  shape  factor;  for  the 

cylindrical  pore  of  circular  cross-section  s*  =  1 .  The  complex  compressibility  is 


C(u»  = 


7  -  (7  - 1) 


Pjj  +  aF  (yfPrH )  /  (iw  Pr) 


(2.15) 


In  this  case,  the  asymptotic  expansions  of  p{u)  an<X  C  (cu)  are  expressed  as: 


pCu)  ^(xjiuj  +  (l  +  sl/Z)p^, 

(2.16) 

C  (a;)  =  (1  /  7PJ  {7  -  (7  - 1)  [wPo  Pr/  a]} 

(2.17) 
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-  for  low  frequencies. 


l  +  s^[a/2iup^) 


1/2 


(2.18) 


(7(0;)  « 


l  +  s^(7-l)(a/2io;Prpo)^^" 


(2.19) 


-  for  high  frequencies. 

Here  the  relation  between  the  structural  factor  and  the  structural  constant  k^,  which  is 
measured  in  experiments,  has  the  simple  form  s^  =  {2/k^^  . 


With  the  approach  it  is  possible  to  determine  the  model  parameters  and  establish  the  relation  of 
the  porous  coating  impedance  at  high  (ultrasonic)  frequency  with  the  impedance  of  the  same 
coating  of  larger  thickness  at  low  (sonic)  frequencies. 

In  general,  the  wave  resistance  and  the  propagation  constant  depend  on  the  scale  parameter  H  as 


fAH) 


fAH)} 


1/2 


(2.20) 


c 


(2.21) 


where  (if)  and  ^  (H)  are  dimensionless  functions  of  the  complex  density  and  compressibility; 
i.e.,  W'  —  const ,  and  r  ~  <a  at  H=const.  For  the  input  impedance  of  the  porous  layer 

H-const  H=const 

attached  to  a  rigid  wall,  we  obtain 


z 

1  '  1 

f  ^  A 

—  ^ 

-St 

II 

r,—h 

pc 

1 

\H=const 

\  J 

,  at  /9W/r7=  , 


(2.22) 


where  and  To  are  acoustic  parameters  at  p^. 

Thus,  the  specific  input  impedance  of  the  porous  layer  of  thickness  h  at  ,Pp  equals  to  the  input 

UJ 

impedance  of  the  porous  layer  of  thickness  h  =  h-^  at  O),  p.  The  latter  can  be  measured  in 

(jj 

experiments.  This  similarity  is  valid  at  /f=const,  that  is  equivalent  to  the  condition  pu  =  p^uj^ 
(assuming  that  the  viscosity  is  constant).  The  wind-tunnel  conditions  (p^  -  O.Olp^^^, 
~  SOOikFfz  ,h~\  mm)  can  be  simulated  in  the  benchmark  experiment  with  p  ~  0.3p  atm  n 

uj  ~  lOkHzm  using  the  porous  layer  of  h'  ~  30/i . 
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2.2.  Measurement  of  propagation  constant,  acoustic  resistance,  and 
sound-absorption  coefficient 

The  propagation  constant  r„  =  +  ia^  and  the  dimensionless  wave  resistance 

W  =  W  -\-iW  are  measured  using  the  method  of  acoustic  interferometer.  This  method  allows 

us  to  obtain  the  absorption  coefficient  for  acoustic  waves  normally  incident  on  the  surface,  and  to 
find  the  complex  impedance.  However,  due  to  geometrical  and  dissipation  limits,  this  method 
can  be  applied  to  disturbances  of  the  sound  frequency  band  only. 


Fig.  n.4  Scheme  of  acoustical  interferometer. 


The  interferometer  (Fig.  n.4)  comprises  a  cylindrical  tube  with  rigid  walls.  One  end  of  the  tube  is 
a  rigid  massive  wall  (tube  plug).  At  the  other  end,  there  is  a  loud-speaker  generating  acoustic 
waves.  A  sample  of  porous  material  is  placed  near  the  tube  plug.  Inside  the  interferometer,  there 
is  a  probe  of  4  mm  diameter  with  a  microphone.  This  probe  can  be  moved  along  the  tube. 
Maximums  and  minimums  of  acoustic  pressure  PmaxtFmin  measured 

along  the  tube  axis.  Theses  data  are  used  for  calculations  of  Fm  and  .  Parameters  of  the 
interferometer  are  given  in  the  table: 


Tube  diameter,  mm 

21 

Effective  frequency  band,  kHz 

0.6 -8.5 

Effective  range  of  pressure,  torr 

760  n  30 

Sound  intensity,  dB 

Probe  phase  correction,  rad 

Probe/tube  cross-section  ratio 

0.03 
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The  first  acoustic-pressure  maximum  is  observed  at  the  surface  of  the  tube  plug,  the  first 
minimum  n  at  the  distance  X./4  from  the  wall,  where  X  is  acoustic  wavelength.  These  properties 
are  used  for  tuning  of  the  interferometer. 


If  the  porous  sample  is  placed  near  the  tube  plug,  the  first  acoustic-pressure  minimum  is  located 
at  a  distance  differing  from  X/4.  The  phase  shift  6  of  the  reflection-from-sample  factor  is 
calculated  using  the  measured  values  of  (at  x^.^)  (Rzhevkin,  1960;  Rosin,  1972; 


Bogolepov,  1986)  as  25  = 


A/2 


,  and  the  standing-wave  factor  is  n  = 


Using  the  values  of  5  and  n,  the  input  impedance  of  the  porous  layer  is  determined  as 

Z^=X^+iY^,  (2.23) 

[n^  -l)sin2^ 


= 


2n 


(n^ +l)-(n^ -l)cos2(^’  “  +l)-(«^ -l)cos2<5 


(2.24) 


The  input  impedance  is  expressed  in  terms  of:  the  porous  sample  thickness  s,  the  back-boundary 
impedance  and  W 


Z  jv  +WshT„s 

“  Z^shT^s  +WchT„s  ’ 


(2.25) 


Since  the  tube  plug  is  absolutely  rigid,  the  back-boundaiy  impedance  is  Z^=oo.  Then,  the  input 
impedance  of  the  porous  layer  is  =  WcthT^s .  If  the  porous  layer  is  placed  at  the  distance 

X/4  from  the  tube-plug  surface,  then  =  WthT^s . 


The  input  impedances  Z„^,^  and  are  measured  for  the  porous  sample  of  the  thickness 

5  being  attached  to  the  rigid  tube  plug  and  being  shifted  to  the  distance  A/4,  respectively.  Then, 
the  acoustic  resistance  is  calculated  as  W  =  ;  its  dimensionless  quantity  is 

W  =-^  =  W  +iW  .  In  turn,  the  propagation  constant  is  determined  as 

l^ffi  ^^^^\l^ex(s+X/4)  ^  ^ex{s)  ' 

s 

It  is  also  possible  to  obtain  and  W  by  measuring  the  input  impedances  and 

the  two  samples  with  the  thickness  s  and  2s  respectively  (both  samples  are  attached  to  the  tube 
plug).  In  this  case,  the  acoustic  resistance  and  the  propagation  constant  are  calculated  using  the 
relations 


W  = 


7 

ctW„s 

m 


= — arch 
2s 


^ex(s)  ^ 
^ex(2s)  ) 


(2.26) 


We  can  also  determine  the  absorption  coefficient  a  for  acoustic  waves  of  normally  incidence, 
using  the  standing-wave  factor  n  and  the  relation 
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(2.27) 


4n 

a  =  — -  . 

+2n  +  \ 


The  upper  limit  of  acoustic  frequency,  ,  corresponds  to  the  case  when  the  half  wavelength  is 

approximately  equal  to  the  largest  transverse  size  of  the  tube;  where  d  is  tube 

diameter.  Below  this  frequency,  the  tube  can  be  treated  as  a  long  line  and  acoustic  waves  are 
almost  flat.  The  lower  limit  is  estimated  from  the  condition  that  the  standing  wave  has  at  least 

c 

two  minimums  and  two  maximums  of  the  acoustic  pressure:  f,^  =  — ^ ,  where  L  is  the 


2L 


interferometer  tube  length. 


To  determine  the  active  component  of  the  material  impedance,  and  its  dependence  on  the  mean 
pressure,  the  resonance  method  is  used.  In  this  case,  the  porous  material  sample  is  placed  at  a 
certain  distance  from  the  rigid  tube  plug.  The  volume  of  air  in  the  tube  between  the  sample  and 
the  wall  is  the  resonator  volume.  The  resonant  frequency  is  observed  when  the  elastic  and  inertial 
components  of  the  system  impedance  are  equal  to  each  other.  At  this  frequency  the  absorption  is 
maximal. 


If  a  quarter  of  the  wavelength  equals  the  sample-to-rigid  wall  distance,  the  air  impedance  in  the 
resonatoris  volume  is  zero;  i.e.,  the  total  wave  load  is  applied  to  the  sample.  In  this  case,  the  real 
parts  of  reactive  (R)  and  inertial  {M)  components  of  the  impedance  are  related  to  the  porous 
material.  If  the  half  wavelength  equals  the  sample-to-rigid  wall  distance,  then  the  absorption  is 
minimal. 


In  Fig.  n.5,  the  absorption  coefficient  of  the  perforated  sheet  is  shown  as  a  function  of  the 
dimensionless  frequency  parameter  kl,  where  k  is  wavenumber  and  /  is  the  distance  to  the  rigid 
wall. 


Fig.  n.5  The  absorption  coefficient  of  the  perforated  sheet  as  a  function  of  the  dimensionless 
frequency  parameter  kl  {k  is  wavenumber,  /  is  the  distance  to  the  rigid  wall);  black  symbols  fi 
/=24  mm;  white  symbols  n  /=50  mm. 
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The  measurements  were  performed  at  atmospheric  pressure.  This  distribution  indicates  the  active 
component  of  the  specific  impedance  is  close  to  1  at  the  resonant  fi'equency;  i.e.,  full  absorption 
is  observed.  The  reactive  component  of  the  impedance  is  not  zero,  since  the  resonant  frequency  is 
lower  than  kill. 

Figure  n.6  shows  the  real  and  imaginary  parts  of  the  specific  input  impedance  for  the  porous 
sample  attached  to  the  wall.  Figure  n.7  shows  the  adsorption  coefficient  a  as  a  function  of  sound 
frequency.  Note  that  below  the  frequency  of  2  kHz  the  absorption  was  too  low  to  be  measured  by 
the  interferometer. 


Fig.  n.6  The  specific  input  impedance  for  the  porous  sample  attached  to  the  tube  plug. 


Fig.  n.7  The  adsorption  coefficient  for  the  porous  sample  attached  to  the  tube  plug. 


The  absorption  coefficient  increases  almost  linearly  in  the  sonic  frequency  band.  The  absorption 
occurs  with  a  small  negative  phase  decreasing  with  frequency.  The  behavior  of  the  input 
impedance  indicates  that  the  elastic  components  (which  are  inversely  proportional  to  the 
frequency)  are  dominant. 

For  the  porous  sample  placed  at  the  distance  A  /  4  from  the  tube-plug  surface,  the  impedances  of 
the  air  layer  is  negligible  compared  to  the  impedance  of  the  rigid  tube  plug.  With  the  appreciable 
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absorption  in  the  volume  between  the  sample  and  the  wall  this  is  not  entirely  correct 
(Voschukova,  1985).  Nevertheless,  this  approach  is  commonly  accepted.  For  this  configuration, 
the  input  impedance  is  the  impedance  of  load  on  the  porous  material.  The  results  are  shown  in 
Figs,  n.8  and  11.9. 


Fig.  n.8  The  specific  input  impedance  for  the  porous  sample  at  the  distance  A/4  from  the  tube 
plug. 


f.kHz 


Fig.  n.9  The  absorption  coefficient  for  the  porous  sample  at  the  distance  A/ 4  from  the  tube- 
plug  surface. 


Theoretically,  as  the  frequency  decreases  the  active  component  of  the  input  impedance  tends  to 
the  specific  mechanical  resistance  and  the  active  component  tends  to  ik\copoh.  The  linear 
approximation  of  the  experimental  data  gives  R  =  0.55,  which  is  very  close  to  the  specific 
mechanical  resistance=0.56,  which  was  measured  directly.  The  parameter  k\  is  1.363,  which  is 
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close  to  the  theoretical  quantity  4/3  for  the  cylindrical  pore  of  circular  cross-section.  Moreover, 
the  effective  pore  diameter,  calculated  from  the  experimental  data,  is  53  pm,  which  is  close  to  the 
average  diameter  50  pm  of  actual  pores. 


Sample 

Structural  factor  s 

Tortuosity 

Added  mass 
coefficient  k\ 

Effective  pore 
diameter,  pm. 

Perforated  sheet 

1.04 

1 

1.363 

53 

At  high  frequencies,  the  active  component  of  the  input  impedance  and  the  absorption  coefficient 
are  significantly  higher  than  their  theoretical  values.  Voschukova  (1985)  pointed  out  that  the 
considerable  increase  of  the  active  component  is  associated  with  the  increase  of  absorption  in  the 
resonator  volume  with  decreasing  of  the  resonator  length. 

The  experimental  values  of  the  input  impedance  allow  us  to  obtain  the  frequency  dependencies 
for  the  wave  resistance  and  the  propagation  constant.  The  latter  can  be  used  for  predictions  of  the 
input  impedances  of  various  absorption  systems  with  this  porous  material.  For  the  sonic 
frequency  band,  the  specific  values  of  the  wave  resistance  and  the  product  Th  of  the  propagation 
constant  and  the  sample  thickness  are  shown  in  Figs.  n.lO,  n.11.  The  solid  and  dotted  lines 
represent  the  theoretical  dependencies,  which  agree  well  with  the  experimental  data  (symbols)  in 
the  sonic  frequency  band  at  atmospheric  pressure.  Using  the  theoretical  model  one  can 
extrapolate  these  dependencies  to  other  frequencies  and  ambient  pressures. 


Fig.  n.lO  The  specific  wave  resistance  of  the  porous  material;  symbols  fi  experiment,  lines  - 
theory. 
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0  2  4  6 

f.kHz 

Fig.  n.12  The  input  impedance  at  atmospheric  pressure;  the  porous  sample  is  at  the  distance  /=24 
mm  from  the  tube-plug  surface;  symbols  -  experiment,  lines  -  theory. 

To  determine  the  dependence  of  the  acoustic  absorption  on  ambient  pressure,  measurements 
were  performed  at  low  air  densities.  The  pressure  was  varied  in  the  range  760-180  mm  m.c. 
Figure  11.12  shows  the  input  impedance  at  atmospheric  pressure,  with  the  porous  material  being 
placed  at  the  distance  /=24  mm  from  the  end  wall.  In  Fig.  11.13,  the  absorption  coefficient  is 
shown  as  a  function  of  frequency  at  various  pressures.  The  frequency  band  covers  the  region,  in 
which  the  resonator  impedance  is  dominated  by  elastic  forces,  as  well  as  the  resonant-frequency 
vicinity.  The  experimental  data  (symbols)  are  compared  with  the  theory  (lines).  Within  the 
considered  pressure  range,  the  active  component  of  specific  impedance  grows  as  the  pressure 
decreases.  This  leads  to  flattening  of  the  resonant  curve.  The  absorption  coefficient  near  the 
resonant  frequency  has  non-monotonic  behavior  versus  pressure.  The  absorption  phase  near  the 
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resonant  frequency  has  a  jump.  As  the  pressure  decreases  this  jump  decreases  and  disappears  at 
the  active  impedance  component=l. 

For  the  sample  locus  /=50  mm,  data  measured  near  the  resonant  frequency  at  various  pressures, 
are  presented  in  Fig.  n.14  (for  the  absorption  coefficient)  and  15  (for  the  phase  shift  of  the 
absorption  coefficient).  Figure  n.  1 6  shows  the  active  component  of  the  input  impedance.  As  in 
the  case  of /=24  mm  (Fig.  n.13),  the  data  (symbols)  agree  well  with  the  theoretical  predictions 
(solid  and  dotted  curves).  This  leads  to  the  conclusion  that  the  theoretical  model  can  be  used  for 
calculations  of  the  porous  layer  characteristics  for  high  frequencies. 


f,kHz 


Fig.  n.13  The  absorption  coefficient  as  a  function  of  frequency  at  various  pressures;  the  porous 
sample  is  at  the  distance  /=24  mm  from  the  tube-plug  surface;  symbols  -  experiment,  lines  - 
theory. 


Fig.  n.14  The  absorption  coefficient  as  a  function  of  pressures  at/=1.24  kHz;  the  porous  sample 
is  at  the  distance  /=50  mm  from  the  tube-plug  surface;  symbols  -  experiment,  lines  -  theory. 
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Fig.  n.15  The  absorption-coefficient  phase  shift  as  a  function  of  pressures  at/=1.24  kHz;  the 
porous  sample  is  at  the  distance  1=50  mm  from  the  tube-plug  surface;  symbols  -  experiment, 
lines  -  theory. 
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Fig.  n.16  The  active  component  of  the  input  impedance  at  various  pressures;  the  porous  sample 
is  at  the  distance  /=50  mm  from  the  tube-plug  surface;  symbols  -  experiment,  lines  -  theory. 
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2.3.  Summary 

The  acoustic  interferometer  has  been  developed  for  benchmark  measurements  of  acoustic 
characteristics  of  porous  materials  within  the  range  of  frequencies  of  0.6  n  8.5  kHz  and  air 
pressures  of  760  fi  30  torr.  In  this  range  of  parameters,  the  wave  resistance  and  propagation 
constant  were  measured  for  the  UAC  of  regular  microstructure,  namely,  the  metal  sheet 
perforated  with  equally  spaced  cylindrical  blind  micro-holes. 

The  theoretical  model,  which  is  based  on  the  imit  problem  describing  propagation  of  acoustic 
disturbances  in  a  long  cylindrical  tube,  was  used  for  prediction  of  the  UAC  impedance  and 
absorption  coefficient.  The  benchmark  measurements  agree  well  with  theoretical  predictions. 

The  theoretical  model  gives  a  similarity  law,  which  allows  for  extrapolation  of  benchmark  data 
to  low  pressures  and  high  frequencies  typical  for  hypersonic  wind  tunnel  tests  and  flight. 

Using  a  combination  of  theory  and  benchmark  experiment,  it  is  feasible  to  design  ultrasonically 
absorptive  coatings  and  evaluate  their  performance  for  hypersonic  laminar  flow  control. 


56 


Chapter  III.  Wind  Tunnel  Experiments  and  Comparison  with 
Theory 

3.1 .  Experimental  model 

The  model  is  a  7  half-angle  sharp  cone  of  500  mm  length.  Its  general  view  is  presented  in  Fig. 
in.2.  Schematics  of  the  model  are  given  in  Fig.  in.2.  The  model  consists  of  the  following  main 
parts: 

•  Sharp  cone  nose  of  65  mm  length 

•  Middle  part  of  65  mm  length  with  the  electric  glow  discharge  source 

•  Base  part  of  370  mm  length  covered  by  the  porous  coating,  which  leading  edge  is  located 
at  1 82  mm  distance  from  the  cone  tip 

•  Unit  for  model  turning 


3.1.1.  Perforated  sheet 


Parameters  of  the  porous  coating  and  manufacturing  are  similar  to  the  GALCIT  T-5  experiments 
(Rasheed  et  al,  2001).  The  laser  drilled  sheet  (see  Fig.  III.3)  of  304SS  stainless  steel  with 
thickness  450  p.m  is  manufactured  by  ActionLaser  Pty  Ltd  (Unit  1,  39  King  Rd,  Hornsby,  NSW 
2077,  Australia).  The  diameter  of  round  cylindrical  microholes  is  50  pm,  the  distance  between 
holes  (spacing)  is  100  ps,  and  the  open  area  (porosity)  is  0.2.  The  total  perforated  area  is 
approximately  33500  mm^  and  the  number  of  holes  is  about  3.35  millions.  Additionally,  the 
rectangular  metal  sheet  (dimensions  109x300  mm)  of  the  same  thickness  and  material  was  made 
solid  (without  holes).  A  40x40  mm^  coupon  of  the  laser-drilled  sheet  was  also  manufactured  for 
benchmark  measurements  of  sound  absorption  (see  Chapter  II). 

The  holes  have  a  conical  shape.  Fig.  in.4,  ni.5  show  two  sides  of  the  perforated  sheet.  The 
average  hole  diameter  is  64±6  pm  at  the  laser  beam  input  (back  side)  and  50±6  pm  at  the  beam 
output  (face  side),  so  that  the  hole  taper  angle  is  about  0.9oo 

The  micrograph  (Fig.  in.6)  shows  that  hole  edges  are  ragged.  The  holes  are  situated  in  rows  with 
spacing  about  100±4  pm.  The  average  distance  between  rows  is  100  pm  (Fig.  HI. 7). 

i 

In  some  areas  of  the  perforated  surface,  the  hole  rows  are  very  close  to  each  other  so  that 
neighboring  holes  are  "merged"  (Fig.  HI. 7).  Four  "merged"  rows  were  found  on  the  total  coating. 
The  mechanical  strength  of  the  perforated  sheet  perpendicularly  to  these  rows  is  low.  The  tensile 
load  is  applied  in  parallel  to  the  "merged"  rows  direction,  so  these  rows  weakly  affected  the 
tension. 
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Fig.  III.3  Perforated  sheet. 
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Fig.  ni.6  Micrograph  of  holes  in  the  perforated  sheet. 


Fig.  in.  7  Perforated  sheet  with  i  merged!  rows  of  holes. 


3.1.2.  Technology  of  manufacturing 

A  cross-section  of  the  model  base  part  is  schematically  shown  in  Fig.  in.8.  Manufacturing  of  this 
part  includes  the  following  steps: 

1 .  Manufacturing  of  a  hollow  metal  truncated  cone. 

2.  Milling  of  through  grooves  of  10  mm  width,  276  mm  length  to  tense  the  perforated  sheet. 

3.  Lathe  turning  of  0.45  mm  depth  hollow  for  the  perforated  sheet. 
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4.  Filing  the  hollow  with  epoxy  onto  the  model  opposite  side. 

5.  Leveling  of  the  epoxy  surface.  To  avoid  misalignment  between  the  cone  surface  and  the 
porous  coating,  the  cone  surface  was  turned  to  the  depth  equal  to  the  coating  thickness. 

6.  Riveting  of  cover  straps  onto  the  sides  of  the  perforated  sheet  for  tension  it  onto  the 
model. 

7.  Cleaning  of  the  perforated  sheet  in  an  ultrasonic  bath  with  tetrachlorethilene. 

8.  Tension  of  the  perforated  sheet  onto  the  model  by  custom  built  tension  mechanism  (Fig. 
III.9).  During  tension  process  a  bend  region  was  close  to  the  peroration  and  solid  strips  at 
the  sides  of  perforated  sheet  were  break  away.  Instead  of  the  broken  strips,  metal  strips  of 
0.25  mm  thickness  were  attached  by  laser  welding.  Because  of  relatively  high  flexibility 
of  the  welded  strips,  the  perforated  sheet  was  successfully  tensed  onto  the  model  (Fig. 

in.io). 

9.  Checking  of  the  porous  coating  level. 

10.  Puttying  and  treatment  of  the  model  surface. 

1 1 .  Final  treatment  of  the  model  with  painting  and  polishing. 

The  roughness  of  the  uncovered  and  polished  part  of  the  cone  is  <0.5  pm. 
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Fig.  EH.  10  Base  part  with  the  regular  porous  coating. 
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3.2.  Experimental  equipment  and  techniques 


3.2.1.  Wind  tunnel 

Experiments  were  carried  out  in  the  hypersonic  blow-down  wind  tunnel  T-326  of  IT  AM  SB  RAS 
(see  Fig.  HI.  12).  This  tunnel  was  designed  for  the  Mach  numbers  Af_=6-14  and  the  unit 
Reynolds  numbers  Rei  =(5•10^-^70T0^)m■*;  there  is  a  set  of  changeable  profiled  and  conic 
axis-symmetric  nozzles.  The  tunnel  is  equipped  with  a  one-step  perforated  ejector. 

The  maximum  pressure  in  the  settling  chamber  is  of  20  MPa;  the  pressure  in  the  ejector  is  of 
0.8  MPa;  the  runtime  is  up  to  30  min.  For  the  Mach  numbers  range  6-10,  an  ohmic  heater  is  used 
for  maintaining  air  temperature  in  the  settling  chamber  up  to  800  K.  The  open-jet  test  section  has 
the  length  of  400  mm  with  the  nozzle  end  diameter  of  200  mm.  To  increase  the  sensor  lifetime 
and  to  exclude  effects  of  dust  particles  on  measurements,  a  block  of  fine  air  cleaning  is  installed 
between  the  air  heater  and  the  nozzle.  This  block  filters  dust  particles  of  the  size  >2  pm.  A 
mechanism  of  model  injection  into  a  free  stream  allows  for  protection  of  models  and  sensors 
from  non-stationary  loadings  during  tunnel  starting.  The  test  section  is  equipped  with  two  optic 
windows  of  220  mm  clear  diameter,  which  are  used  for  model  examination  and  flow 
visualization. 

Experiments  were  performed  at  the  freestream  Mach  number  M«>=  5.95.  The  Mach  number  flow- 
field  nonuniformity  is  0.7%  in  the  flow  core.  The  test  core  diameter  is  approximately  1 80  mm. 
The  noise  level  is  about  1%  that  is  typical  for  conventional  hypersonic  wind  tunnels.  The 
stagnation  pressure  Po  and  stagnation  temperature  To  were  measured  during  the  run  and  kept  at 
constant  with  accuracy  0.06%  and  0.25%,  respectively.  The  pressure  Pq  was  measured  with  the 
help  of  the  pressure  gage  IPD-89008  (the  measurement  ranges  are  1.0  MPa,  1.6  Mpa  and 
2.5  MPa).  A  nameplate  error  of  the  pressure  gage  IPD-89008  is  0.06%  of  the  measurement  range. 
All  measurements  of  Po  were  carried  out  within  the  range  of  1 .6  MPa.  The  stagnation 
temperature  To  was  measured  by  the  k-type  thermocouple  with  the  sensitivity  coefficient  of 
4.024 10’^  mV/grad.  The  accuracy  of  7b  measurement  was  0.1%. 

3.2.2.  Alignment  of  the  model  at  zero  angles 

The  alignment  of  the  cone  along  the  flow  axis  was  achieved  using  the  following  procedure.  To  be 
adjusted  first  was  the  pylon  of  the  model.  The  rotation  unit  with  mirror  mounted  on  its  face  was 
installed  on  the  pylon  (see  Fig.  HI.  13).  The  inclination  angle  of  the  mirror  to  the  mounting-face 
plane  of  the  mirror-rotating  unit  could  be  finely  adjusted.  Further,  the  settling  chamber  of  the 
wind  tunnel  should  be  removed.  Two  transparent  circular  diaphragms  with  1.5  mm  diameter 
central  orifice  were  inserted  into  the  inlet  and  outlet  sections  of  the  nozzle.  A  laser  beam  was 
directed  through  these  openings.  Displacing  the  pylon  frame  we  could  hit  the  mirror  center  by  the 
laser  beam.  Then,  the  mirror  unit  is  rotated  with  simultaneous  adjustment  of  the  inclination  angle 
of  the  mirror  to  the  mounting-face  plane  of  the  mirror-rotating  unit  until  the  beam  reflected  from 
the  mirror  no  longer  wobble.  After  that,  the  mirror  plane  was  assumed  to  be  normal  to  the  axis  of 
rotation  and,  hence,  to  the  axis  of  the  drive  opening  of  the  pylon.  Then,  the  diaphragm  installed 
at  the  nozzle  exit  was  removed.  Varying  the  pylon  installation  angle  we  forced  the  reflected  beam 
fall  onto  the  laser  output  window.  Now,  the  axis  of  the  rotation  unit  could  be  believed  coincide  to 
the  nozzle  axis,  and  the  pylon  was  to  be  rigidly  fixed  to  the  pylon  frame,  as  well  as  the  pylon 
frame  to  the  base  of  the  test  section  of  the  wind  tunnel.  Subsequently,  the  model  was  installed  on 
the  rotation  unit,  and  the  settling  chamber  was  mounted  in  its  proper  place.  This  alignment 
procedure  allows  us  to  install  the  model  at  zero  angle  of  attack  with  accuracy  0.05“ 
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3.2.3.  Traverse  system 

The  T-326  wind  tunnel  is  equipped  with  a  three-component  traverse  system,  which  provides 
0.01  mm  accuracy  in  motion  over  the  coordinates  X,  Y,  Z  (the  coordinates  X,  Y  are  shown  in  Fig. 
m.ll).  The  model  is  turned  around  the  longitudinal  axis  using  the  rotation  unit  providing  0.1° 
accuracy  of  the  turning  angle. 

3.2.4.  Hot-wire  anemometer  system 

Boundary-layer  disturbances  are  measured  using  the  constant  current  anemometer  (CCA) 
manufactured  in  IT  AM.  The  CCA  consists  of  the  two  main  circuits; 

•  Low-frequency  (Ofil  0  Hz)  circuit  to  provide  constant  mean  resistance  of  the  hot-wire. 

•  High-frequency  (<500  kHz)  compensating  amplifier  to  overcome  degradation  of  the  hot¬ 
wire  sensitivity  with  increasing  of  frequency. 

The  measurements  are  carried  out  with  a  single-component  hot-wire  sensor  made  of  tungsten 
wire  of  1  mm  length  and  5  ^im  diameter. 
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Fig.  in.12  Scheme  of  the  T-326  hypersonic  wind  tunnel  of  ITAM  SB  RAS. 

1-  ejector,  2-  diffuser,  3-model,  4-  open-jet  test  section,  5-optical  windows,  6-  axisymmetric 
nozzle,  7-  prechamber,  8-heater,  9, 11 -control  valves,  10-device  to  put  model  into  free  stream. 
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Fig.  in.  13  Alignment  at  zero  angles  procedure. 


3.2.5.  System  of  measurement  control  and  data  acquisition 

The  system  of  measurement  control  and  data  acquisition  structurally  consists  of  the  two  parts: 
1)  System  of  free-stream  parameters  control,  mean  parameters  (P,  TQ  )  measurement  and 
traverse  system  device  operating;  2)  Computer  hot-wire  measurements  system. 

1)  System  of  freestream  parameters  control,  mean  parameters  measurement  and  traverse 
system  device  operating.  The  freestream  parameters  are  determined  from  measurements  of  Po 
and  Tq.  To  measure  Po,  To,  and  Pw,  the  ADC  PCI-91 14HG  (©  ADLINK  Technology  Inc., 
Taiwan)  PC  data  acquisition  card  is  used.  This  ADC  provides  the  following  advanced  features: 

•  16  differential  analog  inputs; 

•  16  bits  high  resolution  AD  conversion; 

•  Sampling  rate  up  to  1 00  kHz; 

•  Input  range:  ±10V,  ±1  V,  ±0.1  V,  ±0.01  V. 

Signals  from  each  sensors  and  thermocouples  are  digitized  with  the  frequency  of  30  kHz.  Each 
set  of  512  samples  are  averaged  and  mean  values  are  stored  in  the  personal  computer  memory. 

The  traverse  gears  are  controlled  using  a  personal  computer  or  manually.  The  software  developed 
for  this  part  of  the  measurement  system  allows  an  experimentalist  to  monitor  free-stream 
parameters  during  the  run,  operate  the  traverse  gears,  record  the  experimental  data  and  start  up 
the  hot-wire  system. 

2)  Computer  hot-wire  measurements  system  consists  of  a  personal  computer  and 
PCI-9812  (©  ADLINK  Technology  Inc.,  Taiwan)  PC  data  acquisition  card.  This  ADC  provides 
the  following  advanced  features: 

•  4  single-ended  analog  input  channels; 

•  12-bit  analog  input  resolution; 

•  Up  to  20MHz  A/D  sampling  rates; 
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•  4  A/D  converters,  one  converter  for  each  analog  channel; 

•  Trigger  Sources:  Software,  Analog  threshold  comparator  using  internal  D/A  to  set  trigger 
level,  and  External  digital  trigger. 

During  hot-wire  measurements  of  natural  disturbances,  1  time-series  of  401408  samples  each  are 
acquired  at  each  measurement  station.  During  measurements  of  artificial  disturbances,  98  time- 
series  of  4096  samples  each  acquired  at  each  measurement  station.  Analog  signal  of  main 
generator  initiating  the  glow  discharge  is  used  for  triggering  ADC.  Since  the  output  signal  of 
CCA  has  a  high  level  of  high-frequency  signals,  the  sampling  rate  is  increased  up  to  5  MHz  to 
avoid  aliasing. 

The  software  allows  for  processing  data  (spectra,  amplitude  and  phase  characteristics  of  artificial 
disturbances)  in  the  real  time  and  plotting  them  on  the  monitor.  During  the  run,  parameters  of 
artificial  disturbances,  sensor  locations  and  free-stream  characteristics  were  stored  into  a  PC  hard 
disk.  At  the  same  time,  the  experimental  program  is  controlled  by  the  second  computer,  and 
exchange  of  information  between  computers  is  performed  by  Ethernet,  protocol  TCP/IP.  Thus, 
the  experimental  process  is  fully  automated.  All  measured  time-traces  are  stored  in  PC  and 
burned  to  CD-ROMs.  This  allows  for  carrying  out  detailed  analysis  of  the  data  after  experiments. 

3.2.6.  The  system  of  high-frequency  glow  discharge 

The  glow-discharge  system  consists  of  a  clock  generator  and  high-voltage  generator.  The  clock- 
generator  signal  is  used  to  initiate  the  high-voltage  generator  and  synchronize  hot-wire 
measurements  with  the  high-voltage  generator  starting  up.  The  high-voltage  generator  produces 
voltage  pulses  up  to  2000  V,  of  1  psec  duration  and  the  pulse  frequency  up  to  400  kHz.  During 
the  run  the  glow  discharge  is  controlled  visually  through  optic  windows  and  a  mirror  in  the  test 
section. 

Ring  and  point  sources  were  tested  in  hypersonic  flows.  During  testing  the  electrode  shapes  were 
optimized  to  improve  the  discharge  efficiency.  Also  the  new  circuit  (providing  more  power, 
higher  frequency  range  and  less  electric  interference  with  hot-wire)  was  successfully  examined. 

3.2.7.  Processing  of  wave  characteristics  of  the  disturbances 

To  obtain  the  amplitude  A  and  phase  O  of  the  controlled  disturbances,  the  discrete  Fourier 
transformation  is  used 

A(x,r,ff)e^(™  =^f,P  ,  (3.i) 

where  N  is  the  samples  count  in  the  time-series,  a)  is  frequency  of  controlled  disturbances, 
pU„(X,  Y,0,tj)  is  the  time-series  of  mass-flow  pulsations. 

To  compare  the  experimental  and  theoretical  data,  the  artificially  excited  wave  packet  is 
decomposed  to  waves.  For  this  purpose,  the  transversal  wave  spectra  were  computed  as 

SA{x,p)e‘^''^^’^^=  (3.2) 

~0Q 

where  SA,  SFnare  the  amplitude  and  phase  spectra  with  respect  to  the  transversal  wave 
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number  p. 

3.2,8.  Measurements  of  natural  and  artificial  disturbances 

Fig.  in.  14  schematically  shows  the  measurement  procedure,  which  includes  the  following  steps: 

•  Vertical  distribution  measurements:  The  hot  wire  is  moved  up  vertically.  The 
measurements  are  carried  out  with  =0.03  mm  step  at  the  distance  of  0.9  mm,  then  with 
=0.05  mm  step  at  the  distance  of  1  mm,  and  then  with  0.1  mm  step  of  at  the  distance  of 
2  mm. 

•  Transversal  distribution  measurements:  The  hot  wire  is  moved  vertically  down  to  the 
pulsation  maximum.  The  maximum  locus  is  determined  at  the  previous  step.  Then  the 
model  is  rotated  within  the  angle  range  0=  ±(20n30)“  depending  on  wave  packet  width, 
and  measurements  are  carried  out  with  the  step  of  =0.5“  This  procedure  also  allows  for 
checking  of  the  boundary-layer  uniformity  in  the  transverse  direction. 

•  Longitudinal  distribution  measurements:  The  model  is  turned  back  to  the  position 
0=0.  Then,  the  hot  wire  is  moved  upstream  parallel  to  the  cone  surface  (at  7“angle)  with 
the  step  AX=0.25  mm  to  the  distance  of  35  mm.  This  procedure  allows  for  measurements 
of  the  artificial  wave-packet  phase  along  the  cone  surface  in  the  layer  of  maximum  mass- 
flow  pulsations 

•  Then  the  procedure  is  repeated 


Fig.  in.  14  Procedure  for  measurements  of  natural  and  artificial  disturbances. 
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3.3.  Results  of  measurements 

Experiments  are  carried  out  at  the  conditions: 

Po=1000  kPa;  ro=385-400  K;  =5.95;  Reioo=(11.5-12.3)xl0®m'';  7y7’o=0.80-0.84. 


3.3.1.  Mean  flow 

Measurements  of  the  mean  boundary-layer  characteristics  are  carried  out  using  the  hot-wire 
anemometer  described  in  Section  3.2.4.  The  upstream  end  of  the  investigated  region  on  both 
solid  and  porous  sides  is  located  ahead  of  the  UAC  leading  edge,  which  is  Jir=182mm 
(Reex=2.810®). 

The  mean  mass  flow  pU,  mean  flow  velocity  U  and  RMS  mass-flow  pulsations  are  shown  in  Fig. 
in.15  n  Fig.  in.20  as  functions  of  the  vertical  coordinate  for  the  solid  and  porous  sides  at  various 
X-stations.  In  Fig.  in.21  n  Fig.  ni.26,  the  same  distributions  are  presented  in  nondimesional 
forms.  The  mass  flow  and  velocity  are  normalized  using  their  values  at  the  upper  boundary-layer 
edge.  The  vertical  coordinate  y  is  normalized  with  respect  to  the  boundary-layer  thickness  5. 

For  X=  138  n  287  mm,  the  profiles  of  pU/(pU)e,  U/Ue  on  the  solid  and  porous  sides  are  similar  to 
each  other;  they  are  close  to  the  correspondent  profiles  calculated  for  the  laminar  boundary  layer. 
The  profiles  of  RMS  pulsations  are  typical  for  the  hypersonic  laminar  boundary  layer.  They  have 
a  large  narrow  peak  in  the  vicinity  ofjp=0.86.  These  data  confirm  that  boundary  layers  are  laminar 
on  the  solid  and  porous  sides. 

Fig.  in.27  compares  U  and  <pU>lpU  on  the  solid  and  porous  sides.  These  profiles  are  close  to 
each  other;  i.e.,  the  porous  coating  does  not  affect  mean-flow  characteristics  and  RMS  mass-flow 
pulsations. 

3.3.2.  Natural  disturbances 

Spectra  of  mass-flow  disturbances  were  measured  at  the  coordinate  y  relevant  to  maximum 
pulsations.  Hot-wire  measurements  are  conducted  at  equally  spaced  X-stations.  The  first  three 
stations  are  located  upstream  from  the  UAC  leading  edge  (X=182  mm).  The  disturbance  spectra 
are  shown  in  Fig.  ni.28a  and  Fig.  III.28b  for  the  solid  and  porous  sides,  respectively.  To  illustrate 
development  of  natural  disturbances,  their  amplitudes  versus  /  and  Rcex  are  show  by  contour 
(Fig.  ni.31)  and  surface  (Fig.  in.32)  plots. 

At  the  three  upstream  stations,  the  spectra  are  similar  on  the  solid  and  porous  surfaces,  as  shown 
in  Fig.  in.29.  Downstream  from  the  UAC  leading  edge,  the  spectrum  behaviors  are  quite 
different: 

•  On  the  solid  side,  the  disturbance  spectra  look  very  similar  to  those  measured  by  Stetson  & 
Kimmel  (1992a)  on  a  sharp  cone  at  M„=6;  they  correspond  to  the  second-mode  instability.  At 
the  first  upstream  station,  the  second  mode  is  observed  at  the  frequency  /  »  430-450  kHz.  Its 
amplitude  quickly  increases  downstream,  whereas  its  central  frequency  decreases  to  270  kHz 
at  the  last  station.  In  the  middle  of  the  measurement  region  (Reex^3,4210®)  the  second  mode 
disturbances  become  lager  than  the  first-mode  disturbances  associated  with  the  frequency 
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band  50-200  kHz;  i.e.,  the  first  mode  grows  slower  than  the  second  mode.  Measurements  at 
the  last  downstream  station  (on  both  solid  and  porous  sides)  indicate  the  beginning  of 
boundary-layer  turbulization. 

•  On  the  UAC  side,  the  second  mode  grows  quickly  in  the  first  three  stations  (upstream  of  the 
UAC).  Then,  this  growth  is  slowed  down  on  the  porous  coating.  Over  the  entire  measurement 
region,  the  second-mode  amplitudes  are  smaller  than  the  first-mode  ones.  Amplifications  of 
the  first-mode  disturbances  are  about  the  same  on  both  solid  and  porous  surfaces. 

To  determine  amplifications  of  natural  disturbances  by  the  boundary  layer,  the  spectra  measured 
in  the  boundary  layer  should  be  referenced  to  the  spectra  of  external  disturbances  in  free  stream. 
However,  the  level  of  high-frequency  (>100  kHz)  freestream  disturbances  is  very  low  at 
hypersonic  speeds  (lower  than  the  hot-wire  anemometer  noise  level).  Therefore,  measurements  of 
these  disturbances  are  not  feasible.  For  this  reason,  disturbance  spectra  are  referenced  to  the 
initial  spectra  measured  in  the  boundary  layer  at  the  first  station.  The  initial  spectra  are  shown  in 
Fig.  in.30  for  the  solid  (line  1)  and  porous  (line  2)  sides.  Since  these  spectra  are  close  to  each 
other,  we  use  for  scaling  one  initial  spectrum  (line  3),  which  is  determined  by  averaging  of  the 
spectra  1  and  2.  To  diminish  effects  of  random  fluctuations,  the  initial  spectrum  Ao  is  smoothed 
by  polynomial  using  the  least  squares  technique  (line  4). 

The  disturbance  spectra  A/Ao  (normalized  by  the  initial  spectrum)  are  shown  in  Fig.  in.33a  and 
Fig.  in.33b  for  the  solid  and  porous  sides,  respectively.  The  normalized  amplitudes  versus  /and 
Reex  are  also  shown  by  contour  (Fig.  111.34)  and  surface  (Fig.  III.35)  plots.  It  is  clearly  seen  that 
natural  disturbances  amplify  in  a  wide  frequency  band  (30-450  kHz)  on  both  solid  and  porous 
sides: 

•  The  first-mode  disturbances  have  nearly  equal  growth  on  the  both  sides.  Maximum 
amplifications  are  observed  in  the  range  100nl30kHz;  they  are  3.3  and  3.6  for  solid  and 
porous  sides,  respectively. 

•  The  second-mode  disturbances  amplify  essentially  stronger  than  the  first-mode 
disturbances.  Their  maximum  amplifications  are  6.9  and  20.4  on  the  porous  and  solid  sides, 
respectively. 

Comparing  these  data  we  conclude  that  the  UAC  of  regular  structure  stabilizes  the  second 
mode  and  weakly  affects  the  first  mode  under  natural flow  conditions. 

3.3.3.  Artificial  wave  packets 

Artificial  disturbances  were  generated  in  the  boundary  layer  at  the  frequency  275  kHz.  At  this 
frequency  the  second  mode  has  maximum  amplitude  in  the  natural  disturbances  spectra  measured 
on  the  solid  side. 

Fig.  in.36  and  Fig.  in.37  show  transversal  distributions  of  the  amplitude  and  phase  of  wave 
packets.  These  distributions  are  nearly  symmetric.  They  are  similar  in  all  cross  sections.  The 
amplitude  distributions  are  maximal  at  the  wave-packet  center.  The  phase  distributions  are  very 
fiat.  The  phase  difference  between  the  neighboring  stations  is  approximately  80  rad.  Fig.  in.38 
shows  the  amplitude  wave  spectra  with  respect  to  the  transversal  wave  number  These  data  are 
obtained  using  the  discrete  Fourier  transform  of  the  transversal  distributions  shown  in  Fig.  in.36 
and  Fig.  in.37.  The  whole  wave  packet  lies  within  the  range  of  >9=±0.5  rad/degree,  which 
approximately  corresponds  to  the  range  of  the  wave-vector  inclination  angle  of  ±20°.  The 
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amplitude  maximums  are  observed  at >^0.  Dimensionless  phase  speeds  Cx=C\lU\  are  shown 
in  Fig.  in.38.  They  weakly  depend  on  the  longitudinal  coordinate  X  and  practically  coincide  on 
the  porous  and  solid  surfaces.  These  features  indicate  that  the  wave  packets  predominantly 
consist  of  the  second  mode  disturbances. 

Amplitudes  of  natural  and  artificial  disturbances  on  the  solid  and  porous  sides  are  shown  in  Fig. 
in.39.  The  disturbance  frequency  is  275  kHz.  On  the  porous  side,  the  amplitudes  of  artificial  and 
natural  disturbances  practically  coincide.  Therefore  the  natural  disturbances  at  this  frequency  are 
predominantly  two-dimensional  waves  of  the  second  mode.  It  is  cleraly  seen  that  the  UAC  leads 
to  substantial  decreasing  of  the  amplitude  growth:  the  maximum  amplitude  is  approximately  3 
times  lower  than  in  the  case  of  solid  wall. 

3.4.  Comparison  with  linear  stability  theory 

Stability  calculations  are  conducted  for  two-dimensional  disturbances  (y9  =  0)  of  the  second 
mode  at  the  flow  parameters  relevant  to  the  experimental  data  shown  in  Fig.  III.39.  Namely:  the 
local  Mach  number  is  M^  =  5325;  the  stagnation  temperature  is  ro=390K;  the  wall 
temperature  is  =  320  K  that  corresponds  to  the  wall  temperature  ratio  7^  =  5.47 1 ;  the  local 

unit  Reynlds  number  is  Re^,  =15.2x10*  m’’;  the  disturbance  frequency  f  =  215  kHz  that 
corresponds  to  the  frequency  parameter  F  =  1.3926x10“’ .  The  fluid  is  a  perfect  gas  of  7  =  1.4 
and  Pr  =  0.708 .  The  viscosity-temperature  dependency  is  approximated  by  Sutherlandis 
formula  (1.7);  the  second  viscosity  is  /i*  =  0.8/i*.  The  spatial  growth  rate  is  calculated  using  the 

relation  (1.30);  the  nonparallel  effects  are  included  as  described  in  Section  1.1.2. 

For  the  porous  wall,  stability  calculations  are  performed  using  the  boundary  conditions  (1.46)- 
(1.48),  in  which  the  characteristic  impedance  and  the  propagation  constant  are  given  by  (1.49)- 
(1.52);  i.e.,  the  rarefaction  effects  are  taken  into  account.  The  UAC  parameters  are  specified  as: 
the  pore  radius  =25  pm,  the  pore  spacing  s’  =100  pm,  the  porosity  ^  =  0.2,  the  porous- 

layer  thickness  h"  =  450  pm,  the  UAC  leading  edge  X  =  \82  mm. 

Figure  III.40  compares  theoretical  amplification  curves  (solid  lines)  with  experimental  data 
(symbols),  which  are  also  shown  in  Fig.  111.39.  In  these  calculations,  the  amplitudes  are  adjusted 
to  the  experimental  values  at  the  initial  station.  The  theory  agrees  well  with  the  experiment  on 
the  solid  side  (compare  the  black  line  with  black  triangles).  The  amplification  curve  calculated 
for  the  porous  surface  (red  line)  passes  above  the  experimental  points  (red  triangles  and  blue 
squares);  i.e.,  the  theory  underestimates  the  UAC  stabilization  effect.  One  of  possible  reasons  of 
this  discrepancy  is  a  conical  shape  of  actual  pores.  The  green  line  shows  the  amplification  curve 
calculated  for  the  pore  radius  averaged  over  the  pore  length  (rj  =  28.5  pm).  With  this  correction 

the  agreement  of  theory  with  experiment  is  satisfactory.  Another  reason  may  be  nonlinear 
effects,  which  may  slow  down  the  second  mode  amplification.  As  mentioned  in  Section  3.3.2, 
the  nonlinear  breakdown  of  natural  disturbances  seems  to  occur  at  the  last  station.  Since  the 
amplitudes  of  artificially  excited  wave  packets  are  significantly  larger  than  those  of  natural 
disturbances,  the  nonlinear  effect  may  be  exposed  earlier.  This  assumption  is  consistent  with 
decreasing  of  the  wave-packet  amplitude  at  the  last  station  (see  blue  squares  in  Fig.  III.40). 

As  shown  in  Figure  in.41,  the  theoretical  phase  speeds  agree  well  with  the  phase  speeds 
measured  for  artificially  excited  wave  packets  on  both  solid  and  porous  sides. 
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3.5.  Summary 


A  model  of  a  T  half-angle  sharp  cone  with  the  UAC  of  regular  microstructure  was 
manufactured.  The  cone  has  a  longitudinal  half  of  its  surface  solid  and  the  other  a  porous  sheet 
perforated  with  equally  spaced  blind  cylindrical  holes  of  the  average  pore  diameter  50  jim,  depth 
450  |J,m  and  average  spacing  100  pm.  The  cone  was  tested  in  the  T-326  hypersonic  wind  tuimel 
at  freestream  Mach  number  5.95.  Stability  measurements  in  the  boundary  layer  on  the  porous 
and  solid  cone  surfaces  were  conducted  using  hot-wire  anemometry.  Spectra  of  natural  and 
artificially  excited  wave  packets  were  obtained  at  various  cross  sections  of  the  cone.  The  wave 
packets  were  generated  by  periodic  glow  discharge  of  a  fixed  frequency. 

Measurements  of  mean  profiles  and  RMS  mass-flow  pulsations  showed  that  the  boundary  layer 
was  laminar  on  the  solid  and  porous  sides  at  all  cross-sectional  stations.  The  porous  coating  did 
not  affect  the  mean  flow  and  RMS  pulsations;  i.e.,  the  UAC  roughness  did  not  cause  premature 
tripping  of  the  boundary  layer. 

Under  natural  conditions,  the  first  (low  frequency)  and  second  (high  fi-equency)  modes  were 
observed  in  the  boundary  layer.  Analysis  of  the  natural-disturbance  spectra  showed  that  the 
second  mode  was  a  dominant  instability  on  the  solid  side.  The  UAC  stabilizes  the  second  mode 
and  weakly  affects  the  first  mode  that  is  consistent  with  the  theoretical  predictions. 

To  evaluate  the  UAC  stabilization  effect,  the  wave  packet  was  artificially  exited  at  the  fi-equency 
275  kHz  relevant  to  the  second-mode  instability.  It  was  found  that  the  two-dimensional  wave  is 
a  dominant  component  of  the  wave  packet.  The  amplification  of  this  component  practically 
coincides  with  the  amplification  of  natural  disturbance  of  the  same  frequency.  This  confirms  that 
natural  disturbances  of  the  high-fiequency  band  are  predominantly  two-dimensional  waves 
relevant  to  the  second-mode  instability.  It  was  shown  that  the  UAC  leads  to  substantial 
decreasing  of  the  wave-packet  growth:  the  maximum  amplitude  on  the  porous  side  is 
approximately  3  times  lower  than  that  on  the  solid  side. 

Stability  calculations  were  conducted  for  the  second-mode  two-dimensional  wave  of  frequency 
275  kHz  and  compared  with  experimental  data.  The  theoretical  amplification  curve  agrees  well 
with  the  experiment  on  the  solid  side.  The  amplification  curve  calculated  for  the  porous  surface 
passes  above  the  experimental  points;  i.e.,  the  theory  underestimates  the  UAC  stabilization 
effect.  With  the  correction  accounting  for  conical  shape  of  actual  pores,  the  agreement  is 
satisfactoiy.  It  is  assumed  that  the  discrepancy  between  theory  and  experiment  in  the 
downstream  measurement  station  is  due  to  nonlinear  effects. 

In  summary,  the  LST  model  discussed  in  Chapter  II  satisfactorily  predicts  the  UAC  stabilization 
effect  and  can  be  used  for  evaluations  of  the  UAC  performance  in  connection  with  hypersonic 
laminar  flow  control. 
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Symbols  1-4  correspond  to  X=146.68,  195.14,  209.33  and  244.65  mm. 
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Fig.  ni.15  The  mean-flow  velocity  profiles  in  the  boundary  layer  on  the  solid  side. 
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Fig.  in.  16  The  profiles  of  pU  in  the  boundary  layer  on  the  solid  side. 
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Fig.  in.  17  The  RMS  pulsations  profiles  on  the  solid  side. 
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Symbols  1-4  correspond  toX=145.56,  180.68,  215.7  and  251.12  mm. 
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Symbols  1-2  correspond  toX=287.2,  321.75  mm 
Fig.  in.  18  The  mean  flow  velocity  profiles  in  the  boundary  layer  on  the  porous  side. 
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Fig.  in.  19  The  profiles  of  pU'm  the  boundary  layer  on  the  porous  side. 


3 


Symbols  1-4  correspond  to  RCexSO  ^^2.18, 2.91,  3.30  and  3.79 


Symbols  1-3  correspond  to  ReexSO  ^“4.23,  4.62,  4.97 


Fig.  in.21  Normalized  profiles  of  U  on  the  solid  side. 
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Symbols  1-4  correspond  to  Reex20^^2.18,  2.91,  3.30  and  3.79 


Symbols  1-3  correspond  to  ReexS^  ^~4.23,  4.62, 4.97 
Fig.  in.22  Normalized  profiles  ofpC/on  the  solid  side. 
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Symbols  1-4  correspond  to  RcexSO  *“2.18,  2.91,  3.30  and  3.79 


Symbols  1-3  correspond  to  RoexS®  *^4.23,  4.62, 4.97 
Fig.  in.23  Normalized  RMS  pulsations  profiles  on  the  solid  side. 
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Symbols  1-4  correspond  to  RcexS^  ^^2.22, 2.81,  3.35  and  3.91 


Symbols  1-2  correspond  to  RcexS®  *“4.50,  5.05 


Fig.  in.24  Normalized  profiles  of  U  on  the  porous  side. 
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Symbols  1-4  correspond  to  ReexS0'*='2.22,  2.81, 3.35  and  3.91 


Symbols  1-2  correspond  to  ReexS®  *=4.50,  5.05 
Fig.  in.25  Normalized  profiles  of pf/on  the  porous  side. 
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Symbols  1-4  correspond  to  Reex20^“2.22, 2.81, 3.35  and  3.91 


Symbols  1-2  correspond  to  Reex20’^=4.50,  5.05 
Fig.  III.26  Normalized  RMS  pulsations  profiles  on  the  porous  side. 
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Fig.  in.27  Comparison  of  U  and  <pU>lpU  distributions  on  the  solid  side  (1  fi  Reex=2.1820‘^) 
and  on  the  porous  side  (2  n  Reex=2.2220  *). 


84 


Amplitude  Amplitude 


0  100  200  300  400  500 

f,  kHz 


Fig.  in.28  Natural  disturbances  spectra  in  the  boundary  layer  on  the  solid  (a)  and  porous  (b) 
surfaces. 
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Fig.  in.29  Comparison  of  spectra  measured  upstream  from  UAC;  lfl3  fl  solid  side,  Reex  lO' 
=2.28,  2.55,  2.84;  4n6  n  side  with  UAC,  Reex  l0'*=2.22,  2.50,  2.78. 


Fig.  in.30  The  initial  spectra. 

1  fi  solid  side  (Reex  lO"^  =2.28);  2  n  side  with  UAC  (ReexTO'^  =2.22);  3  n  average  spectrum; 
4  -  smoothed  average  spectrum. 


Fig.  in.31  Contour  plot  of  natural-disturbance  spectra  for  solid  (a)  and  porous  (b)  surfaces. 
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Fig.  in.32  Surface  plot  of  natural-disturbance  s 
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Fig.  in.34  Contour  plot  of  normalized  natural-disturbance  spectra  for  solid  (a)  and  porous  (b) 
surfaces. 
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Fig.  III.36  Transversal  distributions  of  the  artificial  wave-packet  amplitude  on  the  porous  side. 
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Fig.  in.37  Transversal  distributions  of  the  artificial  wave-packet  phase  on  the  porous  side. 
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Fig.  III.39  Downstream  development  of  the  second-mode  disturbances  of  frequency  275  kHz. 

1  n  natural  disturbances  on  solid  side;  2  fi  natural  disturbances  on  porous  side;  3  fi  artificial 
disturbance  at  y^O,  porous  side 


Fig.  111.40  Comparison  of  theoretical  amplification  curves  (lines)  with  experiment  (symbols)  for 
porous  and  solid  sides;  disturbances  of  frequency  275  kHz. 
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Conclusions 


Experimental  and  theoretical  studies  of  hypersonic  boundary-layer  stabilization  using  a  passive 

ultrasonically  absorptive  coating  (UAC)  of  regular  microstructure  were  conducted  in  order  to 

evaluate  the  UAC  performance  for  hypersonic  laminar-flow  control.  The  results  of  these  studies 

are  summarized  as  follows: 

1.  Linear  stability  problem  for  three-dimensional  disturbances  in  two-dimensional  boundary 
layers  on  a  porous  wall  was  formulated  in  the  framework  of  the  asymptotic  method  of 
multiple  scales.  This  problem  differs  from  the  typical  stability  problem  because  of  new 
boundary  conditions  associated  with  absorption  of  disturbance  energy  by  a  porous  coating. 
The  boundary  conditions  are  formulated  in  terms  of  the  porous-wall  acoustic  admittances, 
which  depend  on  the  UAC  microstructure. 

2.  Acoustic  properties  of  a  cylindrical  pore  were  analyzed  in  the  framework  of  linear  acoustic 
theory.  Dependencies  of  the  dynamic  density  and  dynamic  compressibility  on  gas  pressure 
and  disturbance  frequency  were  derived  including  the  gas  rarefaction  effects.  The  analysis  is 
focused  on  the  most  practical  case  when  the  molecular  mean  free  path  is  smaller  (or  much 
smaller)  than  the  pore  radius.  Calculations  indicate  that  acoustic  properties  of  the  perforated 
panel  essentially  depend  on  Knudsen  number  in  practical  ranges  of  pressure  and  frequency. 

3.  Parametric  calculations  of  the  UAC  effect  on  the  second-mode  stability  are  performed  for  the 
T-326  wind  tunnel  conditions  including  the  gas  rarefaction  effect.  In  the  case  of  adiabatic 
wall,  the  net  stabilization  effect  is  small  or  even  negative  for  disturbances  of  relatively  small 
frequencies.  The  UAC  stabilization  dramatically  increases  as  the  wall  temperature  ratio 
decreases.  This  trend  is  preferable  for  practice  because  surfaces  of  actual  hypersonic  vehicle 
are  maintained  at  low  wall-temperature  ratios. 

4.  The  rarefaction  effect  is  appreciable  for  the  T-326  conditions.  The  presence  of  the  Knudsen 
layer  on  the  pore  wall  leads  to  deeper  penetration  of  disturbances  into  pores  and  provides 
stronger  stabilization  of  the  second  mode. 

5.  The  acoustic  interferometer  has  been  developed  for  benchmark  measurements  of  acoustic 
characteristics  of  porous  materials  within  the  range  of  frequencies  of  0.6  n  8.5  kHz  and  air 
pressures  of  760  fi  30  torr.  In  this  range  of  parameters,  the  wave  resistance  and  propagation 
constant  were  measured  for  the  UAC  of  regular  microstructure,  namely,  the  metal  sheet 
perforated  with  equally  spaced  cylindrical  blind  micro-holes. 

6.  The  theoretical  model,  which  is  based  on  the  unit  problem  describing  propagation  of  acoustic 
disturbances  in  a  long  cylindrical  tube,  was  used  for  prediction  of  the  UAC  impedance  and 
absorption  coefficient.  The  benchmark  measurements  agree  well  with  theoretical  predictions. 

7.  The  theoretical  model  gives  a  similarity  law,  which  allows  for  extrapolation  of  benchmark 
data  to  low  pressures  and  high  frequencies  typical  for  hypersonic  wind  tunnel  tests  and  flight. 
Using  a  combination  of  theory  and  benchmark  experiment,  it  is  feasible  to  design 
ultrasonically  absorptive  coatings  and  evaluate  their  performance  for  hypersonic  laminar  flow 
control. 

8.  A  model  of  a  7“  half-angle  sharp  cone  with  the  UAC  of  regular  microstructure  was 
manufactured.  The  cone  has  a  longitudinal  half  of  its  surface  solid  and  the  other  a  porous 
sheet  perforated  with  equally  spaced  blind  cylindrical  holes  of  the  average  pore  diameter  50 
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)xm,  depth  450  p,m  and  average  spacing  100  |J,m.  The  cone  was  tested  in  the  T-326  hypersonic 
wind  tunnel  at  freestream  Mach  number  5.95.  Stability  measurements  in  the  boundary  layer 
on  the  porous  and  solid  cone  surfaces  were  conducted  using  hot-wire  anemometry.  Spectra  of 
natural  and  artificially  excited  wave  packets  were  obtained  at  various  cross  sections  of  the 
cone.  The  wave  packets  were  generated  by  periodic  glow  discharge  of  a  fixed  fi'equency. 

9.  Measurements  of  mean  profiles  and  RMS  mass-flow  pulsations  showed  that  the  boundary 
layer  was  laminar  on  the  solid  and  porous  sides  at  alt  cross-sectional  stations.  The  porous 
coating  did  not  affect  the  mean  flow  and  RMS  pulsations;  i.e.,  the  UAC  roughness  did  not 
cause  premature  tripping  of  the  boundary  layer. 

10.  Under  natural  conditions,  the  first  (low  frequency)  and  second  (high  frequency)  modes  were 
observed  in  the  boundary  layer.  Analysis  of  the  natural-disturbance  spectra  showed  that  the 
second  mode  was  a  dominant  instability  on  the  solid  side.  The  UAC  stabilizes  the  second 
mode  and  weakly  affects  the  first  mode  that  is  consistent  with  the  theoretical  predictions. 

11.  To  evaluate  the  UAC  stabilization  effect,  the  wave  packet  was  artificially  exited  at  the 
frequency  275  kHz  relevant  to  the  second-mode  instability.  It  was  found  that  the  two- 
dimensional  wave  is  a  dominant  component  of  the  wave  packet.  The  amplification  of  this 
component  practically  coincides  with  the  amplification  of  natural  disturbance  of  the  same 
frequency.  This  confirms  that  natural  disturbances  of  the  high-fi'equency  band  are 
predominantly  two-dimensional  waves  relevant  to  the  second-mode  instability.  It  was  shown 
that  the  UAC  leads  to  substantial  decreasing  of  the  wave-packet  growth;  the  maximum 
amplitude  on  the  porous  side  is  approximately  3  times  lower  than  that  on  the  solid  side. 

12.  Stability  calculations  were  conducted  for  the  second-mode  two-dimensional  wave  of 
frequency  275  kHz  and  compared  with  experimental  data.  The  theoretical  amplification  curve 
agrees  well  with  the  experiment  on  the  solid  side.  The  amplification  curve  calculated  for  the 
porous  surface  passes  above  the  experimental  points;  i.e.,  the  theory  underestimates  the  UAC 
stabilization  effect.  With  the  correction  accounting  for  conical  shape  of  actual  pores,  the 
agreement  is  satisfactory.  It  is  assumed  that  the  discrepancy  between  theory  and  experiment 
in  the  downstream  measurement  station  is  due  to  nonlinear  effects. 

In  summary,  the  LST  model,  benchmark  measurements  of  acoustic  absorption  by  the  porous 
coating  samples  and  stability  measurements  on  a  sharp  cone  in  Mach=6  wind  tunnel  confirmed 
the  concept  of  hypersonic  boundary-layer  stabilization  using  passive  porous  coatings  of  regular 
microstructure.  Theoretical  and  experimental  results  of  this  effort  can  be  used  for  evaluation  of 
the  UAC  performance  and  design  of  hypersonic  laminar  flow  control  systems. 
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Future  Effort 

The  following  issues  need  to  be  addressed  in  further  studied: 

1.  The  wall  temperatures  of  actual  hypersonic  vehicles  are  substantially  lower  than  the  adiabatic 
wall  temperature.  Stability  calculations  showed  that  the  UAC  performance  dramatically 
increases  as  the  wall  temperature  ratio  decreases.  These  theoretical  predictions  need  to  be 
verified  by  stability  experiments  on  the  model  with  cooled  wall.  Such  experiments  can  be 
conducted  in  the  T-326  Mach=6  wind  tunnel  of  ITAM. 

2.  Roughness  of  the  porous  coating  may  lead  to  premature  tripping  of  the  boundary  layer. 
Detailed  studies  of  this  detrimental  effect  on  stability  and  transition  are  required  to  formulate 
tolerances  on  pore  size  and  spacing. 

3.  It  remains  to  be  shown  that  the  second-mode  stabilization  leads  to  the  increase  of  laminar 
rati.  At  present,  we  have  only  indirect  experimental  evidence  in  favor  of  this  connection. 
Namely,  the  experiments  in  the  GALCIT  T-5  shock  tunnel  of  Caltech  (Rasheed  et  al.,  2001) 
showed  significant  increase  of  the  transition  Reynolds  number  on  the  porous  surface.  Further 
experimental  studies  are  needed  to  clarify  this  issue. 

4.  Systematic  benchmark  measurements  of  acoustic  properties  of  various  porous  materials  may 
help  to  identify  optimal  UAC  microstructures.  Such  measurements  can  be  performed  in 
ITAM  using  the  acoustic  interferometer  discussed  in  Chapter  n  of  this  report. 

5.  Further  refinements  of  the  theoretical  model  are  also  needed  to  address  the  pore  end  effects 
and  treat  acoustic  disturbances  in  porous  materials  of  complex  microstructure.  The  theory 
may  be  verified  by  the  benchmark  experiments  of  Item  4. 
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Appendix 

D  =  dldy,  n'  =  dnldT ,  w  =  2(e-l)/3,  r  =  2(e  +  2)/2,  where  e  =  1.2  corresponds  to  the  ratio 

2 

of  the  second  viscosity  to  the  first  viscosity  //j  ///  =  — e  =  0.8 ;  ;jr  =  [—  +  (aU  -  <»)]’' . 

3  U 


Non-zero  elements  of  the  matrix  Hq  in  (1.14)  are 

Du 


DT  Dfl.  RDU 


=  or"  +  +  KaU  -a)—,  H^=  — ^ =  -ia[(m  + 1)— + — ]  + 


HT 


H^=i—-(rn  +  \)}M^a(aU  -co),  H^={m  +  \)^{aU  -  co)  ■ 

H  T 


T  ^  HT 
D{/l'DU) 

5 

M 


_  /I'DU. 
“26  “  ..  > 


M 


DT  i 

/fj,  =  -ia ,  /fjj  =  — ,  =  -i'iMl  (aU  -d),  H^^=-  (aU  -  d) ,  =  -ip ; 

„  .  .  DT  ^Du.  „ 

^41  =-iOcX(r—  +  2-^),  =-iccx, 

1  fi 

rr  .  2  02  DuDT  D^T  .  R  ,  rr  M 

DT  Du, 

^44  =  -iXriMllocDU  +  (—  +  -^\aU  -  d)] , 

I  n 

H,s=iX[roc^+C(^^  +  r^iaU-d)],  H,,=iXj{aU-d), 

H„=-iPx{.r^  +  2^),  H,,=-iPx-, 

1  // 

DT 

=-2Pr(r-l)M^DU,  =  RPr - 2ia{r-\)M^  PiDU, 

fxr 


=  -iRPx^'^  - d) , 


(aU-d) 


^.^uXDUf  D^n 


Du 


H,,  =a‘^+P^+  iRPx  ^ - ir- \)Ml  Pr ^  =-2^; 

/XT 
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H. 


83 


m 

T 


M 


i^-  (m  +  \)'}MlPiaU  -  CO) , 

II 


H^,={m  +  \)^{oiU-(o),  H^T=a^+p^+i{aU-(0)^,  H^=-^. 

1  fii  /J' 

The  vector  G  =  H,Fo  in  (1 . 1 8)  are  expressed  as 
Fo=(/.,/2v..,/8)", 


G,  =  Gj  =  Gf  =  0 , 

R 


Gi  = 


A^+fo^+ 

T  3xj  T  dy 


MA_A 

T 


A  BT 


U  dT 


T  dx, 


T  dx, 


dx, 


G,  mlf,  -f 


V 


A 


dU  ^  dVp 
dx,  dy 


-jmX 


& 


r  dy  j  T  dy  dy 


f 


V  f 


Ge  = 


i?Pr 


M 


’_A^_ 

fs]( 

T  0^ 

[  T  T^)[ 

dxj  dy  j 


+{r-\)Ml 


dy 


T  dxij 


juT  dy  ’ 
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